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During  the  last  decade,  considerable  progress  has  been  made  in  radar  signal  processing,  and  this 
report  states  its  present  status.  The  three  broad  areas  of  coherent  processing,  noncoherent  detection, 
and  track-while-scan  systems  are  discussed. 

Specifically,  in  the  area  of  coherent  processing  the  subjects  of  sidelobe  cancelers,  adaptive 
radars,  MTIs,  and  doppler  processing  are  discussed.  In  the  adaptive  processing  area,  both  the  maxi- 
mum signal-to-noise  and  least  mean-square  methods  are  described,  and  special  emphasis  is  given  to  — 
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the  problem  of  convergent  rate.  The  moving-target  detector  (MTD)  is  used  as  an  example  of  dopp- 
ler  processing. 


In  the  area  of  noncoherent  detection,  various  integrators  are  discussed.  Among  these  are  the 
moving  window,  feedback  integrator,  two-pole  filter,  binary  integrator,  and  batch  processor. 

Methods  of  obtaining  a constant  false-alarm  rate  using  either  adaptive  thresholding  or  nonparametric 
detectors  are  also  considered. 

A general  outline  of  a track-while-scan  system  is  considered  first.  Then,  detailed  discussions  of 
the  tracking  filter,  maneuver-following  logic,  track  initiation,  and  correlation  logic  are  presented. 
Finally,  methods  of  integrating  data  from  several  radars  are  discussed. 


iTi.sSiJll  ID!  _ 1 

Wktti  fcctU* 

’DC 

l«tt  smmm  □ j 

UHMKOCKCn 

o 

j JDSTIf  15*1101 

||  

IISTMiUTIOMMIUUIUn  **»B 

tut. 

AVAIL  Alt/*  IfBIU  1 

CONTENTS 


INTRODUCTION  I 

COHERENT  PROCESSING  1 

Sidelobe  Cancelers  2 

Adaptive  Arrays  and  Radars \ . 6 

Moving  Target  Indicators  14 

Doppler  Processing  18 

Noncoherent  MTI  Indicators  19 

NONCOHERENT  DETECTION 20 

Classical  Theory 20 

Integrators 22 


Moving  Window  22 

Feedback  Integrator  23 

Two-Pole  Filter  24 

Binary  Integrator  24 

Batch  Processor  24 

False  Alarms  25 


Adaptive  Thresholding 26 

Nonparametric  Detectors 28 


Sequential  Detectors 


30 


TRACKING  SYSTEM  31 

System  Outline  32 

Tracking  Filters  33 

Maneuver-Following  Logic  35 

Track  Initiation  38 

Correlation  Logic  38 

Radar  Integration  41 

CONCLUDING  REMARKS  44 


I . iii 

L_ 


SURVEY  OF  RADAR  SIGNAL  PROCESSING 


INTRODUCTION 

. 

During  the  last  decade  considerable  progress  has  been  made  in  radar  signal  processing. 
This  progress  is  directly  traceable  to  the  lowered  cost  and  increased  speed  of  digital  hardware 
and  computers  and  to  more  sophisticated  techniques  in  adaptive  processing  and  tracking  sys- 
tems. 

This  survey  of  radar  signal  processing  will  neglect  waveform  design  and  include  the 
t track-while-scan  systems.  Waveform  design  will  be  neglected  because  it  has  received  consider- 

able attention  elsewhere,  with  the  books  of  Rihaczek*  and  Cook  and  Bernfeldt  covering  the 
subject  in  detail.  On  the  other  hand,  although  track-while-scan  systems  properly  fall  under  the 
heading  of  radar  data  processing,  it  does  not  make  sense  to  have  an  automatic  detection  system 
unless  it  is  accompanied  by  a tracking  system.  Therefore,  since  tracking  is  a necessary  part  of 
the  entire  system,  the  survey  will  include  it. 

Thus  this  survey  of  radar  signal  processing  will  consider  the  three  broad  areas  of  coherent 
processing  (processing  of  amplitude  and  phase),  noncoherent  processing  (processing  of  ampli- 
tude), and  track-while-scan  systems.  The  subjects  will  be  discussed  in  the  same  order  as  the  radar 
signal  passes  through  the  radar  system.  Specifically,  in  the  area  of  coherent  processing  the  sub- 
jects of  sidelobe  cancelers,  adaptive  antennas,  and  MTIs  (moving-target  indicators)  will  be 
covered.  In  the  area  of  noncoherent  detection,  methods  of  obtaining  a constant  false-alarm 
rate  (CFAR)  using  either  adaptive  thresholding  or  nonparametric  detectors  will  be  emphasized. 
The  section  on  the  tracking  system  will  cover  the  tracking  filter,  correlation  logic,  track  initia- 
tions, maneuver-following  logic,  and  a basic  overview  of  an  entire  tracking  system. 

COHERENT  PROCESSING 

In  the  area  of  coherent  processing,  adaptive  processing  will  receive  considerable  attention. 
There  are  two  approaches  to  adaptive  processing:  the  method  of  maximum  signal-to-noise  ratio 
(S/N)  due  to  Howells*  and  Applebaum§  and  the  least-mean-square  method  (LMS)  due  to 
Widrow  and  Hoff#.  The  two  methods,  although  appearing  quite  different,  yield  almost 
equivalent  results.  So  that  both  methods  will  be  presented,  the  LMS  method  will  be  used  dur- 
ing discussion  of  sidelobe  cancelers,  and  the  method  of  maximum  S/N  will  be  used  during  dis- 
cussion of  adaptive  arrays  and  radars.  For  adaptive  radars  special  consideration  will  be  given  to 


*A  W K ihaczek,  Principles  of  Hmh-Resolulion  Radar.  McGraw-Hill,  New  York.  1969. 

tC\  E.  Cook  and  M.  Bernfeld,  Radar  Signals.  A n Introduction  to  Theory  and  Appluanon.  Academic  Press.  New  York.  1967. 
tP.  W.  Howells,  IEEE  Trans  Anlennas  and  Propagation  AP-24.  575-584  (1976) 

§S.  P.  Applebaum,  IEEE  Trans.  Antennas  and  Propagation  AP-24,  585-598  (1976). 

» B Widrow  and  M.  E Hoff.  IRE  WF.SCON  Conv.  Rec  , 96-104,  I960. 
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the  problem  of  convergent  rate.  Finally,  MTIs  will  be  discussed  and  the  moving-target- 
detector  (MTD)  system  will  be  used  as  an  example  of  doppler  processing. 

Sidelobe  Cancelers 

The  basic  idea  of  a sidelobe  canceler  (a  device  that  attempts  to  eliminate  interference 
entering  through  the  antenna  sidelobes)  is  shown  in  Fig.  1.  The  signal  S of  interest  enters 
through  the  main  lobe  of  the  antenna,  and  the  jamming  (interfering  signal),  which  is  much 
stronger  than  the  signal  of  interest,  enters  through  the  sidelobe  of  the  main  antenna.  The  aux- 
iliary antenna  is  an  omnidirectional  antenna,  and  it  will  be  assumed  that  the  signal  entering 
the  omnidirectional  antenna  is  much  smaller  than  the  jamming  Ja  and  can  be  neglected,  since 
the  signal  and  jamming  now  have  the  same  antenna  gain.  (The  treatment  of  the  signal  in  the 
auxiliary  channel  can  be  found  in  Widrow  et  al.*)  The  adaptive  filter  produces  an  output  F 
which  is  as  close  as  possible  to  the  input  jamming  7.  The  filter  output  is  then  subtracted  from 
the  main  input,  producing  an  output  Z = S + J — Y.  If  the  filter  output  is  an  exact  replica  of 
7,  the  output  is  the  desired  signal  S. 


MAN 


Fig.  1 — Concept 
of  adaptive  noise 
canceling 


The  filter  is  controlled  by  adjusting  its  parameters  to  minimize  the  output  power.  To 
show  that  this  minimization  will  force  Y to  be  a replica  of  7,  a development  in  Widrow  et  al.*  is 
repeated.  First,  assume  S,  7,  and  Ja  are  zero-mean  random  variables,  5 is  uncorrelated  with  7 
and  7a,  and  Ja  (and  hence  F)  is  correlated  with  7.  The  expected  output  power  is 

£|Z2)  = E{S2)  + £1(7  - F)2)  + 2£1S(7  - F)}  = £{S2}  + £|(7  - Y) 2}.  (1) 

Adjusting  the  filter  to  minimize  £(Z2)  is  equivalent  to  minimizing  £|  (7  — F)2),  since  Fis  un- 
correlated with  5;  that  is,  F is  the  best  least-squares  estimate  of  the  jamming  7.  Furthermore, 
since  Z — S = J — Y,  minimizing  £(  (7  — F)2)  causes  Z to  be  the  best  least-squares  estimate 
of  the  signal  5. 

The  adaptive  filter  for  obtaining  a least-squares  estimate  of  a desired  signal  S can  be 
described  by  a weighting  vector  IF,  where  IFr  = (IF,,  IF2 Wn  ) and  Fdenotes  the  tran- 

spose, operating  on  the  input  Ja  = X,  X1  — (.v,,...,.v„ ).  Thus  the  filter  output  is 

F = XTW,  (2) 


'B  Widrow,  J R Glover.  Jr  , i M McCool.  J Kauniu.  C.  S Williams,  R.  H.  Hearn.  J.  R.  Zeidler,  E.  Dong,  Jr.,  and  R 
C Goodlin,  Proc  IF.EE  At.  IA9M7I6  (1975). 
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and  the  error,  defined  as  the  difference  between  the  input  signal  and  the  filter  output,  is 

t = 5 + J - XTW.  (3) 

The  least-mean-square  (LMS)  adaptive  filter  adjusts  the  weighting  vector  W to  minimize  the 
mean-square  error.  The  squared  error  is 

€ 2 = (S  + J)2  - 2(S  + J)XTW  + WTXXTW.  (4) 

Taking  the  expected  value  of  (4),  letting  the  vector  P be  the  crosscorrelation  between  J and  X 
( P = E[JX] ),  and  letting  the  matrix  K be  the  covariance  matrix  of  X (K  = E[XXT J ),  one  ob- 
tains 

f(€2)  - £{52)  + E\J2]  -2 PTW  + WtKW.  , (5) 

To  find  the  minimum  of  (5)  with  respect  to  W,  the  gradient  V of  (5)  is  set  to  zero,  yielding 
the  optimal  weight  vector 

W - K~XP.  (6) 

The  LMS  adaptive  algorithm  is  an  iterative  method  of  finding  an  approximate  solution  to 
(6).  The  algorithm  has  the  advantage  of  not  requiring  an  explicit  measurement  of  the  correla- 
tion function  or  inversion  of  the  covariance  matrix.  Specifically,  the  LMS  algorithm  uses  the 
method  of  steepest  descent  to  solve  (6);  that  is,  the  next  weight  vector  + l is  equal  to  the 
old  weight  vector  plus  a step  in  the  direction  of  the  negative  gradient: 

"0-H  = Wj  -mV,..  (7) 

The  gradient  of  the  squared  error  on  the  y'th  iteration  is 

V,  = Ve,2  = V(5  + J - Xj  Wj  ) 2 = -UjXj.  (8) 

Thus  the  next  weight  is  given  recursively  by 

Wj  + \ = Wj  + I^jXj  (9) 

and  is  known  as  the  Widrow-HofT  LMS  algorithm.  The  parameter  n is  a factor  which  controls 
the  rate  of  convergence  and  the  stability  of  the  method.  It  has  been  shown*t  that  (9)  con- 
verges to  the  optimal  solution  as  long  as  m is  between  zero  and  the  reciprocal  of  the  largest 
eigenvalue  of  the  covariance  matrix  K.  Shown  in  Fig.  2 is  a typical  learning  curve  and  an  aver- 
age of  48  learning  curves  for  the  LMS  algorithm.  The  average  reveals  the  basic  exponential 
nature  of  the  learning  curve.  For  the  radar  case  Xj  represents  the  sample  fromy'th  range  cell; 
consequently  the  number  of  iterations  corresponds  to  the  number  of  range  cells. 

In  principal,  if  the  situation  shown  in  Fig.  1 is  correct  (no  uncorrelated  noise  in  each 
channel  and  no  signal  in  the  auxiliary)  the  jamming  can  be  completely  canceled.  However,  if 
the  situation  is  as  shown  in  Fig.  3,  total  cancellation  cannot  be  accomplished.  Specifically,  the 
performance  of  the  canceler  can  be  described  by  the  ratio  R of  S/N  at  the  output  to  S/N  at  the 


•B  Widrow,  P.  E.  Mantey,  L.  J.  Griffiths,  and  B B.  Goode,  Proc.  IEEE  55,  2U3-2I59  <1967) 
tR  L.  Riegler  and  R T Compton,  Jr.,  Proc.  IEEE  61.  748-758  (1973). 
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Fig.  2 — Typical  learning  curves  for  ihe  LMS  algorithm. 
(From)  B.  Widrow  el  al„  Proc.  IEEE  63.  1692-1716  (1975). 
courtesy  of  the  Institute  of  Electrical  and  Electronics  En- 
gineers.) 


primary  input  (main  antenna).  Widrow  et  al*  have  shown  that  this  ratio  R for  steady  state 
• after  convergence)  can  be  expressed  as 


[A(z)  + 1 ) IB(:)  + 1 ] 
A(z)  + A(z)  B(z)  + fl(z)' 


(10) 


where  A(z)  and  Biz)  are  noise-to-noise  ratios 


A (z)  = Sq  (c)/5„  (z) 

and 


(11) 


fllz)  = 5,  (z )/S„  (z)|  H (z)| 2. 


(12) 


*B  Widrow.  J R Glover.  Jr , J M.  McCool.  J Kaunitz,  C.  S Williams.  R H.  Hearn.  J R Zeidler.  E Dong,  Jr.,  and  R 
C.  Goodlin,  Proc  IEEE  63.  1692-1716  (1975). 


4 


NRL  REPORT  8117 


in  which  S^,  St , and  Sn  are  the  power  density  spectra  of  the  noises  w0,  /ni . and  n respectively 
and  H(z)  is  the  channel  transfer  function  for  the  correlated  noise  (jamming).  It  is  obvious 
from  (10)  that  the  cancellation  is  limited  by  the  uncorrelaled  noise  components  in  the  primary 
and  reference  channels.  When  the  jamming  is  much  stronger  than  the  uncorrelated  noise 
components,  A (r)  and  Biz ) are  small  and 


R Aiz)  + Biz)' 


03) 


giving  a large  improvement  in  the  output  signal-to-jamming  ratio.  However  the  improvement 
indicated  by  (13)  is  rarely  achieved  in  practice.  Factors  limiting  performance  include  the  finite 
time  for  the  adaptive  process,  the  presence  of  signal  components  in  the  auxiliary  channel,  mul- 
tipath problems,  and  misadjustment  caused  by  gradient  estimation  noise  in  the  adaptive  pro- 
cess.* Furthermore,  in  theory  N omnidirectional  antennas  (and  associated  cancellation  loops) 
are  needed  to  cancel  N jammers.  However,  because  of  multipath  propagation,  the  energy  from 
a single  jammer  can  enter  the  antenna  from  several  directions  and  for  all  practical  purposes  ap- 
pears to  be  from  several  jammers.  Therefore  in  practice  one  requires  several  times  as  many 
cancellation  loops  as  jammers. 


Recently  F.  Kretschmer  and  B.  Lewist  have  developed  an  improved  algorithm  for  simula- 
tion of  the  Applebaum-Howells  adaptive  loop  and  for  use  in  adaptive  processing.  The  LMS  al- 
gorithm discussed  above  is  given  by 

Wj  + 1 = Wj  + 2 pLtjXj.  (9) 

This  is  commonly  used  to  simulate  and  analyze  the  Applebaum-Howells  adaptive  loop  in  the 
form 

WJ  + \ — kWj  + G ( 1 -kUjXf.  04) 

where  k = 1 — 1/r,  with  r being  the  filter  smoothing  constant,  and  G being  the  gain  term. 
Thus  in  both  algorithms  the  next  weight  is  derived  in  terms  of  the  present  error  and  sample. 
Kretschmer  and  Lewis  point  out  that  for  fast  loops  + ) as  given  by  (9)  and  (14)  is  not  the 
proper  weight.  Rather,  for  better  cancellation  and  more  realistic  canceler  loop  simulation, 
Wj  + 1 should  be  calculated  from 

ITj  + 1 = Wi  +2mVm*/+i.  <15> 

In  effect,  by  using  the  sample  X:  to  calculate  the  weight  + a phase  shift  is  introduced 
which  can  result  in  loop  instability.  Kretschmer  and  Lewis  have  shown  (for  the  Applebaum- 
Howells  application)  that  the  stability  condition  of  the  LMS  algorithm  is 

| G(\  - k )\Xj  \ 2 - k\  <1  (16) 

and  that  their  improved  algorithm  is  unconditionally  stable. 


Comparison  of  the  LMS  algorithm  with  the  improved  algorithm  was  made  using  comput- 
er simulations.  Correlated  Gaussian  noise  (mean  = 0,  variance  = 2)  was  used  as  an  input  to 
the  main  and  auxiliary  channels  of  the  sidelobe  canceler.  At  the  250th  range  cell  a constant 
signal  at  S/N  = —20  dB  is  introduced.  The  signal  residue  for  both  algorithms  with  canceler 


"8.  Widrow,  P.  E.  Manley,  L J Griffiths,  and  B B Goode.  Proc.  IEEE  55.  2143-2159  (1967). 

♦F  Kretschmer  and  B L.  Lewis,  "An  Improved  Algorithm  for  Adaptive  Processing."  NR1.  Report  8084.  Dec  1976 
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parameters  of  k — 1 2 rr  (0.000124)  and  (7  = 100  is  shown  in  Figs.  4 and  5.  Although  the 

LMS  algorithm  had  unstable  performance,  the  improved  algorithm  had  completely  stable  per- 
formance. Also,  for  slow  loops  there  will  be  ringing  in  the  LMS  algorithm,  which  will  result  in 
degraded  cancellation  performance  In  a previous  paper  Kretschmer*  investigated  cascading 
sidelobe  canceler  stages  as  a method  of  obtaining  improved  cancellation  ratios  and  transient 
responses.  Thus  a higher  effective  loop  gain  would  be  achieved  with  low  actual  loop  gains, 
which  are  required  for  stable  operation.  In  lieu  of  their  later  work,  the  improved  algorithm 
provides  another  way  of  obtaining  high  loop  gains.  Lewis  and  Kretschmer  are  now  working 
on  a open-loop  digital  implementation  of  a sidelobe  canceler. 


Fig.  4 — Adaplive-canceler  response  of 

the  LMS  algorithm 


The  sidelobe  canceler  removes  the  jamming  signal  after  it  has  entered  the  main  antenna. 
Adaptive  arrays,  which  require  individual  receiving  elements,  attempt  to  prevent  jamming  from 
entering  the  antenna  receive  pattern  by  placing  a receiving  antenna  null  in  the  direction  of  the 
jammer.  Before  commencing  with  a discussion  of  adaptive  arrays  and  radars,  it  is  pointed  out 
that  the  September  1976  issue  of  the  IEEE  Transactions  on  Antennas  and  Propagation  is  a spe- 
cial issue  on  adaptive  arrays  and  contains  many  interesting  articles. 

Adaptive  Arrays  and  Radars 

Qualitatively,  in  an  adaptive  array  the  received  signal  is  the  weighted  sum  of  the  signal  at 
the  individual  receiving  elements,  with  the  weights  being  a function  of  the  received  signal. 
The  theory  of  adaptive  arrays  was  first  discussed  by  Applebaum.t  and  VVidrow  et  al.t  have 


•F.  F.  Kretschmer,  IEEE  International  Radar  Conf..  181-185,  1975. 

+S.  P Applebaum,  "Adaptive  arrays,"  Syracuse  University  Research  Corp.  Report  SPL-769.  June  1964. 
»B.  W id  row.  P,  E Manley.  L.  J.  Griffiths,  and  B B Goode,  Proc.  IEEE  55,  2147-2159  (1967). 
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made  major  contributions  to  the  theory,  however  a later  development  of  Brennan  and  Reed* 
will  be  followed.  Their  approach  is  similar  to  Applebaum’s  in  that  they  maximize  S/N,  which 
they  show  is  equivalent  to  maximizing  the  probability  of  detection  when  the  noise  is  Gaussian 
distributed. 

Let  the  radar  be  composed  of  N receiving  elements,  and  let  the  last  M time  samples  from 
each  element  be  processed.  Thus  there  are  n — NM  space-time  samples.  Define  S to  be  a 
complex  (amplitude  and  phase)  //-vector  which  contains  the  desired  signal  components,  and 
define  A"  to  be  a complex  //-vector  containing  the  noise  samples.  The  radar  return  2 is  given 

by 

2 = S + X.  (17) 


To  detect  the  signal  S,  the  radar  output  is  passed  through  a linear  filter  described  by  a weight- 
ing vector  W.  Thus  the  output  of  the  detector  (the  filter)  is 

Y = WTZ.  (18) 


Brennan  and  Reed  showed  that  S/N  at  the  output  of  the  filter  is 


| WTS\  2 
wt'kw' 


(19) 


where  the  asterisk  indicates  the  complex  conjugate  and  K is  the  noise  covariance  matrix, 
K =*  E{X  X r 1,  X having  zero  mean.  Consequently  what  is  required  is  the  value  of  (V  that  max- 
imizes (19).  If  the  Schwarz  inequality  is  used,  it  can  be  shown  that  the  maximum  value  of 
U9)  is  SrE~iS*  and  that  this  value  is  obtained  when 


W = a K~'S*. 


(20) 


'L.  E.  Brennan  .mil  I S.  Heed.  IHEF.  Trans  Aerospace  and  Electronic  Systems  AES-9,  237-252  (1973) 
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where  a'  is  an  arbitrary  nonzero  complex  number.  This  criterion  has  been  known  for  some 
time.*  However,  it  is  rarely  used,  since  K is  not  known  a priori;  and  if  K is  estimated,  it  has 
been  extremely  difficult  to  invert  K in  real  time. 

What  makes  the  Brennan-and-Reed  approach  different  from  other  adaptive  array  process- 
ing is  not  the  ability  to  place  spatial  nulls  in  the  direction  of  jammers  but  rather  the  temporal 
processing  that  is  equivalent  to  a motion-compensated  MTI  (moving-target  indicator).  The 
compensated  MTI  behavior  is  obtained  by  selecting  the  proper  steering  signal  S.  The  selection 
of  the  steering  signal  S will  be  illustrated  for  the  case  of  an  airborne  coherent  pulsed  radar. 


Assume  that  the  return  is  range  gated,  there  are  range  celis,  and  the  return  from  the 
;th  cell  is 


z (y)  = x(j)  + s(j). 


(21) 


The  return  signal  from  the  rth  receiving  element  and  /nth  time  sample  can  be  written  as 


S,(/w)  = bre""v.  r = 1 N. 


(22) 


where  y = ~A-nVTI\  is  the  doppler  phase  shift,  with  V being  the  relative  velocity  of  the  tar- 
get, T being  the  time  between  transmitted  pulses,  and  X being  the  radar  wavelength.  The 
quantity  br  is 


b , = A,e 


i<b,  +iR 


= 1, 


jV, 


(23) 


where  Ar  is  the  signal  amplitude  at  the  rth  element,  8 is  a constant  phase  fac'or,  and  d>r  is  the 
relative  phase  between  the  target  and  the  rth  element.  For  a linear  array  with  element  spacing 
cl,  the  phase  angles  6r  for  a signal  arriving  at  an  angle  <l>  with  respect  to  the  array  normal  are 

<t>r  = — sin  t/»,  r = 1 N.  (24) 

Thus  the  expected  signal  for  a linear  array  can  be  obtained  by  substituting  (23)  and  (24)  into 

(22). 


Both  clutter  and  target  will  have  returns  of  the  form  of  (22).  Since  the  velocity  of  the  tar- 
get (and  consequently  the  relative  velocity  F)  is  unknown,  it  is  impossible  to  specify  Sfor  the 
optimal  weighting  given  by  (20).  However,  since  (22)  is  computable  for  ground  clutter  as  a 
function  of  the  radar-clutter-cell  geometry,  one  selects  a steering  signal  5 which  is  orthogonal  to 
the  ground-clutter  vector  5'.  Thus  the  purpose  of  5 is  to  reject  the  clutter,  not  to  detect  the 
target.  This  is  about  as  close  to  an  optimal  detector  as  one  can  obtain,  since  it  can  be  shown* 
that  no  uniform  most-powerful  test  exists  when  the  target  velocity  is  unknown. 

As  an  example  let  M = 2 and  assume  on.e  wants  to  detect  a target  in  a direction  normal 
to  the  direction  of  the  platform  velocity  (the  radar  is  sidelooking).  Then  Sr(m ) = Are'n , and 
for  uniform  amplitude  taper  (Ar  = 1,  r = 1 N ) the  clutter  signal  is 

S T = e*  [1 1,1 1 ].  (25) 

The  appropriate  steering  signal  S which  is  orthogonal  to  S',  STS  * = 0,  is 

ST  = (I I,  -1 -1  ].  (26) 


*11  L.  VanTrees,  IEEE  Trans  Military  Electronics  MIL-9,  216-229  (1965) 
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which  corresponds  to  a target  at  1/2  the  blind  speed  of  the  radar: 

(2?) 

Thus,  if  (26)  is  used  in  (20),  the  detector  is  optimized  for  canceling  main-beam  clutter.  We 
now  consider  how  (20)  can  be  implemented  adaptively. 

Brennan  and  Reed  use  tne  method  of  steepest  ascent  to  maximize  S/N: 


W'  KW 

The  recursive  algorithm  for  steepest  ascent  is 

W(j  + \ ) « W(j)  + y/z  (J)VF[W(J)).  (29) 

where  V F[W(j)\  is  the  complex  gradient  of  devaluated  at  W(j),  which  has  been  shown  to 
be 


WTS 

Vf  =2  - - S* 

WT  KW 


Wr‘S* 

WrKW 


If  K is  assumed  known  and  m (j)  is  chosen  to  be  a constant,  one  can  apply  known  theorems* 
to  show  W(j)  approaches  a critical  point  as  a limit.  Thus,  if  WTO)  is  sufficiently  close  to  the 
optimal  value,  W(j)  approaches  a'K^S*  in  the  limit. 

The  trouble  with  using  (30)  in  (29)  is  that  Vdis  a nonlinear  function  of  W(j),  which  in 
some  adaptive  systems  can  cause  computational  difficulties.  Hence  the  algorithm  was  linear- 
ized by  noting 

WTS  1 A 

lim  — — = — 7 = a.  (31) 

j-oo  yyT  KW  a 

Thus,  if  n (J)  equals  a constant  /a,  (29)  reduces  to 

W(j  + \)  = W(j)  +/xa[S*  - a*K(j)W(j)\.  (32) 

where  K(j ) is  a statistical  estimate  of  the  unknown  covariance  matrix  K.  The  best 
(maximum-likelihood)  estimate  of  K is 

K(j)  = Z*(j)ZT(j).  (33) 

Brennan  and  Reed  then  showed^  that  (32)  converged.  Specifically,  the  expected  value  of  (32) 
converges  to  a K -1S*,  where  K = E[K(j)\  for  all  j,  if  ZJ_j)  are  independent  and  0 < /a  < 
max  A,,  where  A,  ( / = 1 , ...,  «)  are  the  eigenvalues  of  K. 

The  block  diagram  of  the  adaptive  radar  is  shown  in  Fig.  6,  and  the  implementation  of  an 
adaptive  loop  is  shown  in  Fig.  7.  The  steady-state  antenna  pattern  can  be  calculated  from  (20), 
and  the  S/N  improvement  can  be  found  from  S'K~^S'.  However  in  many  radar  environ- 
ments the  clutter  has  a temporal  and  spatial  variation;  consequently  ihc  rate  of  covergence  is 
important.  To  study  this  phenomena,  computer  simulations  were  used. 

*M.  J D Powell.  SIAM  Rev.  12.  79-97  (1970) 
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The  basic  parameters  for  a ten-element  adaptive  array  using  only  one  time  sample  {N  = 
10  and  XI  — 1)  are  given  in  Table  1.  In  the  first  simulation,  30  discrete  clutter  points  were  uni- 
formly distributed  in  the  two  symmetrical  intervals  [17°,  90°]  and  [ —17°,  — 90°],  and  the  radar 
was  looking  normal  to  the  aircraft  velocity  vector.  The  simulation  results  are  summarized  in  Fig.  8, 
where  the  base  of  the  plot  is  45  dB  below  the  peak  gain.  The  back  antenna  pattern  is  the  ini- 
tial receiving  pattern,  the  middle  eight  patterns  are  from  range  cells  200  to  1600  in  200  range- 
cell intervals,  and  the  last  pattern  is  the  steady-state  pattern.  Since  there  are  30  interference 
sources  and  only  10  elements,  it  is  impossible  to  put  a null  at  each  intereference  angle.  Rather 
the  adaptive  array  follows  two  strategies:  it  widens  the  main  beam  and  consequently  lowers 
the  general  sidelobe  level,  and  it  places  receiver  nulls  at  transmitter  maximums  and  vice  versa. 
After  1600  interactions  all  but  1.6  dB  (27.3  — 25.7)  of  the  maximum  signal-to-clutter  improve- 
ment has  been  obtained. 


In  the  second  simulation  the  30  clu'ter  points  were  placed  nonsymmetrically  about  zero 
in  the  interval  [15°,  45°].  The  simulation  results  are  summarized  in  Fig.  9.  Although  the 
sidelobes  are  reduced  in  the  proper  angular  interval,  after  1600  iterations  only  24.7  dB  of  the 
possible  44.1-dB  improvement  in  the  signal-to-clutter  ratio  has  been  obtained.  Brennan  and 
Reed  have  shown  that  the  time  behavior  of  the  weights  is  a sum  of  exponentials  of  the  form 


W. 


= £ c„ e 

i-i 


l-(CXj-H  )f/r  I 


(34) 


where  r is  the  time  constant  and  G is  the  gain  of  the  low-pass  filter.  Thus  the  rate  of  conver- 
gence is  controlled  by  the  smallest  eigenvalue  of  K\  specifically,  the  effective  time  constant  is 
r/(6'*mjn  + l ).  This  suggests  that  rapid  convergence  can  be  obtained  by  selecting  G to  be 
large  and/or  r to  be  small.  However  this  is  not  a useful  solution  to  the  convergence  problem 
since  Brennan  et  al.*  have  shown  that  the  total  output  noise  power  in  the  adaptive  array  is 


P 


WTKW  1 


G A v 

2t  L 1 
i-i 


(35) 


*L  E Brennan.  E.  1.  Pugh,  and  I S Reed.  IEEE  Trans  Aerospace  and  Electronic  S> stems  AES-7,  254-2P7  '1971) 
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Table  I — Parameters  Assumed  in  a 
Simulation  of  an  Adaptive  Receiving  Array 

Ten-element  linear  array 

Element  patterns  isotropic  over  -jt/2  > v ^ ir/2 
Half-wave-spaced  elements 
Uniformly  illuminated  transmit  array 

30  scatterers  in  the  sidelobe  region,  equally  spaced  in  angle 
No  interference  for  -H  ( < 0 < H 
Each  receiving-element  weight  controlled  adaptively 
Simulation  of  1600  independent  sets  of  input  signals  (range 
resolution  cells) 

No  receiver  noise 


Fig.  8 — Projeciograph  plot  of  ihe  gain 
of  a ten-element  adaptive  array  in  the 
case  of  symmetric  clutter  distribution. 
The  improvement  in  the  signal-to- 
sidelohe  clutter  ratio  from  the  initial  re- 
ceiving pattern  (at  the  read  is  27 .5  dtt 
for  steady  state  (pattern  at  the  front) 
and  25.7  dB  after  1600  iterations. 
(From  L.  E.  Brennan  and  L S.  Reed, 
IEEE  Trans.  Aerospace  and  Electronic 
Systems  AES-9,  237-252  « 1973).  cour- 
tesy of  the  Institute  of  Electrical  and 
Electronics  Engineers.) 


Fig  9 — Projeciograph  plot  of  Ihe  gain 
of  a ten-element  adaptable  array  in  the 
case  of  nonsymmetric  clutter  distribu- 
tion The  improvement  in  the  signal- 
to-sidelobe  clutter  ratio  is  44  1 dB  for 
steady  state  (not  shown)  and  24.7  dB 
after  1600  iterations  (pattern  at  the 
front).  (From  L F.  Brennan  and  L S. 
Reed,  IEEE  Trans.  Aerospace  and 
Electronic  Systems  AF.S-9.  2.77-252 

(1977),  courtesy  of  the  Institute  of 
Electrical  Engineers.) 
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where  W is  the  average  weight  vector  in  the  absence  of  loop  noise  (departure  from  steady 
slate).  The  quantity  \VTAfW  is  the  noise  power  when  W = K~XS*.  Consequently . the  output 
power  has  been  increased  by  the  factor  G Ia,/2t  due  to  loop  noise.  Thus,  when  K contains 
both  small  and  large  eigenvalues,  it  is  impossible  to  select  a G and  r which  yield  both  rapid 
convergence  and  low  loop  noise.  To  avoid  the  convergence  problem.  Reed  el  al.*  have  sug- 
gested a direct  computation  of  the  weights. 

The  maximum-likelihood  estimate  of  K,  assuming  the  noise  is  Gaussian  distributed,  is 


k =y  £ Z*(j)ZT(j). 

L l-\ 


Since  Z*(j)ZT (j)  is  an  «-by-/t  matrix  of  rank  1,  L must  be  for  the  inverse  to  exist.  Then 
the  filter  has  the  form 

W = K~XS\  (37) 

The  output  S/N  for  (37)  normalized  by  the  maximum  S/N,  SrK  -IS*.  which  corresponds  to 
(20),  is 

(STk~'S*)2 

[>(K)  = ; — ; : ; . (38) 

(SrK~xS*)  (SrK~x  KK  “ S*) 

The  expected  value  of  (38)  is 

F\,,(k)\  = (L  + 2 -«)/(£  + !).  (39) 

Thus  the  average  loss  can  be  kept  less  than  3 dB  (£(p  (k))  ^ 1/2)  by  letting  L ^ 2n. 

However,  whereas  the  adaptive  loops  of  Fig.  6 require  //  complex  multiplications,  the 
sample-matrix  inverse  method  requires  approximately  o3  complex  multiplications.  To  reduce 
the  complexity  of  the  method,  one  can  update  the  covariance  matrix  using 

K,  = (1  — a)  Kj  _j  + ocZ*(j)ZT(j).  (40) 

where  n is  the  weight  applied  to  the  current  sample.  Then  the  inverse  of  Kj,  given  Kj  , isf 

• , A/-1,  „ | ['V/-! z-0>]  [zrO'ix,-l| 

K -\  _ — 2 ; ' . (41) 

J I “«  1 “«  (1  -«)  +aZT{j)K~_\  Z*(j ) 

This  method  of  updating  the  inverse  requires  approximately  2 n2  complex  multiplications.  The 
average  computation  time  for  updating  the  weights  W'  depends  on  how  frequently  they  must 
be  updated.  For  example,  depending  on  the  radar  environment,  updating  the  weights  every 
PRF  using  (36)  may  be  quite  adequate;  consequently  the  computation  time  may  be  less  than 
that  of  the  adaptive  loops. 

Brennan  et  al.t  compared  the  convergent  rates  of  the  three  methods  using  a computer 
simulation  illustrating  airborne  MTI  performance.  The  results  of  the  simulation  are  shown  in 


’I  S Recil,  J 1).  Mullctt.  and  t.  E.  Brennan.  IEEE  Trans.  Aerospace  and  Electronics  Systems  AES-10.  S5J-8t>)  (1974) 
tj  M Shapard.  I)  Edelblute.  and  li  Kinnison,  Naval  Undersea  Research  and  Development  Center  Report  NI  C 
I N - 528.  May  1971 

1 1.  I Brennan . J I)  Mallet t.  and  I S Reed.  IEEE  Trans.  Antennas  and  Propagation  AES-24.  607-615  1 19761. 
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F ig.  10.  In  both  instances,  (a)  forward  looking  and  (b)  sidelooking,  the  two  methods  of  calcu- 
lating K 1 provide  an  excellent  convergent  rate.  Figure  10  indicates  an  MTI  gain  of  plus  100 
dB,  but  in  practice  the  MTI  gain  would  be  limited  to  a lower  figure  by  internal  clutter  motion. 

Most  work  on  adaptive  arrays  and  radars  has  been  limited  to  theoretical  studies.  However 
there  has  been  some  experimental  work  at  Ohio  Slate  University,*  the  Naval  Research  Labora- 
tory,! and  the  Wide- Aperture  HF  Radio  Research  Facility  operated  by  Stanford  Research  Insti- 
tute^ 


Moving-Target  Indicators 

Moving-Target  Indicators  (MTIs)  were  first  investigated  in  the  1940's,  and  they  have 
been  discussed  in  detail  in  the  books  by  Skolnik#**  and  Naihansontt-  The  coherent  MTI, 
the  most  common  MTI,  uses  an  internal  coherent  reference  source  to  distinguish  a moving  tar- 
get from  fixed  clutter  returns.  The  MTI  signal  is  obtained  by  coherently  subtracting  the  re- 
turned voltages  from  successive  transmitted  pulses: 

Z,  ( j)  = Z,  (./)  - Z,  (42) 

where  Z,  (j)  is  the  dh  returned  pulse  in  ihc./ih  range  cell.  Larger  clutter  attenuations  can  be 
obtained  by  using  multiple  pulses.  The  frequency  (doppler)  response  of  the  MTI  is  that  of  a 
bandpass  filter. 

The  most  serious  problems  associated  with  MTI  are  limiting  and  blind  speeds.  The  first 
of  these  can  be  covered  very  simply.  In  the  classic  paper  or  Ward  and  Shradertt  it  was  shown 
that  MTI  improvement  could  be  degraded  by  20  dB  in  a three-pulse  canceler  by  limiting  the 
clutter  return.  Their  work  showed  that  the  degradation  was  fundamental  to  limiting  and  that 
consequently  a large  dynamic  range  is  required  to  avoid  limiting. 

The  major  problem  with  MTI  is  that  blind  speeds,  corresponding  to  doppler  frequencies 
higher  than  Nyguist  rate,  occur  at 

= Jf-  1 = '-2-3-  <43) 

Thus  for  an  L-band  (1.3-GHz)  radar  with  a PRF  ot  300  pps  the  blind  speeds  occur  at  multiples 
of  approximately  70  knots.  Because  of  the  width  of  the  clutter  notch  (rejection  region  of  the 
canceler),  many  air  targets  would  not  be  detected.  There  are  several  solutions  to  the  problem 


•R  T.  Compton.  IEEE  Trans  Antennas  anil  Propagation  A P- 2-1.  697- 706  (1976). 

tW.  F.  Gabriel,  "Proceedings  Adaptive  Antenna  Systems  Workshop  March  11-13.  Vol  I"  NRL  Report  7801  Sect 
1974. 

»!.  I Griffiths,  IEEE  Trans  Antennas  and  Propagation  A P-24,  707-720  (1976). 

ST.  W.  Washburn  and  L E.  Sweeney.  Jr.,  IEEF.  Trans  Antennas  and  Propagation  AP-24.  721-7.1’  (1976) 

0 M I.  Skolnik.  lnimtttiction)n  Radar  Systems.  McGraw-Hill,  New  York,  1962 
"M  I Skolnik . editor  Riulur  HonilbmiK.  McGraw-Hill.  New  York.  1970. 
ttf  F..  Nathanson,  Rmlar  lies ■::»  Principles,  McGraw  Hill,  New  Y ork.  1969 
OH  K Ward  and  W W Shrader.  FASCON  Convention  Record.  168-173  1968 
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of  blind  speeds  in  MTIs.  Among  ihese  arc  variable  PRF,  siaggered-PRF  MTI,  and  dual- 
frequency MTI. 

The  simplest  solution  is  to  use  a variable-PRF  system.  If  an  interpulse  period  of  T is 
used,  a blind  speed  of  VH  is  obtained.  Then,  if  the  interpulse  period  is  changed  by  a small 
fraction  /-,  the  blind  speed  changes  by  the  same  fraction  r;  and  the  smallest  common  blind 
speed  is  Ftf/(  1 — r).  Thus,  if  an  L-band  radar  has  two  PRFs.  300  pps  and  270  pps,  the  blind 
speed  of  the  radar  system  is  approximately  700  knots.  There  are  two  disadvantages  of  such  a 
system:  (a)  second-time-around  clutter  (clutter  beyond  the  unambiguous  range,  caused  by- 
ducting  at  sea  or  high-altitude  long-range  clutter  such  as  mountains  and  chaff)  passes  through 
the  MTI,  and  (b)  the  constant  PRF  for  a two-  or  three-pulse  burst  makes  the  system  more 
vulnerable  to  jamming.  The  simple  solution  to  (a),  using  an  extra  filler  pulse  (transmitting 
three  pulses  but  only  using  the  last  pulse  out  of  a two-pulse  MTI),  makes  situation  (b)  worse. 

An  elegant  solution  to  the  blind-speed  problem  is  the  staggered-PRF  MTI.  The  basic 
MTI  configuration  is  shown  in  Fig.  11.  The  interpulse  durations  t,  are  constrained  by  the  rela- 
tion 


-i. 


(44) 


where  FH  is  the  first  blind  doppler  frequency  and  £,  are  integers  for  all  /.  Capon*  showed  that 
the  optimal  weights  |o,l  for  minimizing  the  output  clutter  residue  while  retaining  some  fraction 
of  the  average  gain  of  the  filter  (this  constraint  avoids  the  trivial  solution  a,  =0,  for  all  /)  are 
the  components  of  the  eigenvector  associated  with  the  smallest  eigenvalue  of  the  clutter  co- 
variance  matrix.  This  procedure  ignores  what  happens  in  the  filter  passband.  Hsiao  and 
Kretschmert  developed  a procedure  for  setting  the  interpulse  periods  to  minimize  the  RMS 
passband  ripple  while  maintaining  the  minimum  clutter  residue.  A typical  response  is  shown 
in  Fig.  12.  The  basic  trouble  with  this  system  is  that  second-time-around  clutter  will  not  be 
canceled. 

A third  solution  to  the  blind-speed  problem  is  the  dual-frequency  MTI  first  discussed  by 
Kroszczy nskit§  and  later  by  Hsiao#.  The  system  works  by  transmitting  two  frequencies 
whose  ratio  / is  slightly  less  than  I,  filtering  out  the  sum  signal  and  retaining  the  difference  sig- 
nal. The  system  performance  is  basically  that  of  a low-frequency  radar;  hence  the  blind-speed 
problem  is  reduced.  The  detrimental  factor  is  that  the  clutter  improvement  factor  is  reduced 
by  several  dB.  A typical  filter  response  for  a dual-frequency  MTI  is  shown  in  Fig.  13. 
Although  the  passband  response  is  quite  variable,  no  attempt  has  been  made  to  reduce  the 
variation  by  changing  r.  Hsiao  indicates  that  the  staggered-PRF  MTI  is  preferable  to  the  dual- 
frequency MTI.  However  this  author  believes  that  the  dual-frequency  MTI  should  not  be  dis- 
carded that  readily.  An  alternate  solution,  and  possibly  a better  one,  is  to  operate  individual 
MTIs  at  the  two  frequencies. 


’J.  Capon.  It.lU  Trans  Information  Theory  IT-10.  152-159  (1964). 

tJ  K Hsian  and  F t Kretschmer.  Jr.,  The  Radio  ami  Electronic  Engineer  43.  65  ' 93  (1973). 
tj  Kros/czynski.  Radio  and  Electronic  Engineer  34.  157-159  (1967). 
jj  Kros/c/ynski,  Radio  and  Electronic  Engineer  39,  172-176  (1970) 
tt  I K Ilsiao,  The  Radio  and  Electronic  Engineer  45.  351-356  1 1975) 
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Tig  II  — A staggered-PRF  MTI  filter. 
(Prom  J K Hsiao  and  F.  F.  Kretsch- 
mer. Radio  and  Electronic  Engineer  43, 
689-693  (1973),  courtesy  of  the  Institu- 
tion of  Electronic  and  Radio  En- 
gineers.) 


t 


DOPPLER  FRECKJENCr/PRF 


Fig.  12  — Frequency  response  for  a 

seven-pulse  staggered-PRF  MTI  filter. 
(From  Hsiao  and  Kretschmer,  Radio 
and  Electronic  Engineer  43,  689-693 
(1973).  courtesy  of  the  Institution  of 
Electronic  and  Radio  Engineers  ) 


Fig  13  — Target-signal  gain  function 
of  a dual-frequency  MTI  system  with  r 
= 0.89  (From  1 K Hsiao.  Radio  and 
Electronic  Engineer  45.  351-356  (1975), 
courtesy  of  the  Institution  of  Electronic 
and  Radio  Engineers  ) 
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Doppler  Processing 

An  MTI  canceler  provides  near  opiimal  target  detection  in  clutter  but  provides  little  or  no 
improvement  against  receiver  noise.  McAulay*  formulated  the  problem  as  a classical  detection 
problem  and  showed  that  the  opiimal  detector  could  be  structured  approximately  as  an  MTI 
canceler  followed  by  a narrow-band  doppler  filter  bank.  This  structure  has  the  practical  advan- 
tage of  greatly  reducing  the  dynamic  range  required  at  the  input  of  the  filter  bank.  In  this 
configuration,  the  MTI  canceler  provides  improvement  against  clutter,  and  the  doppler  filter 
bank  provides  improvement  against  noise. 

The  moving-target  detector  (MTD),  developed  by  Lincoln  Laboratory^  for  the  FAA, 
uses  this  type  of  processing.  During  1976  the  MTD  was  tested  with  a modified  FPS-18  radar  at 
the  FAA  facility  in  Atlantic  City,  N.J.  The  modified  FPS-18  radar  is  an  S-band  radar  instru- 
mented to  48  n.mi.  The  range  cell  is  approximately  1/16  n.mi.,  the  beamwidth  is  1.5°,  the  scan 
rate  is  15  rpm,  and  20  pulses  are  returned  as  the  radar  sweeps  past  the  target. 

A block  diagram  of  the  MTD  signal  processor  is  shown  in  Fig.  14.  An  azimuth  cell  is 
defined  as  a half  beamwidth  (0.75°)  and  contains  ten  pulses,  with  the  time  lapse  for  the  ten 
pulses  being  referred  to  as  a coherent  processing  interval  (CPI).  In  a CPI  the  ten  pulses  are 
passed  through  a three-pulse  MTI  canceler,  and  the  eight  output  pulses  (two  pulses  are  needed 
to  load  the  MTI)  serve  as  an  input  to  an  eight-point  FFT,  the  points  being  weighted  to  provide 
low  filter  sidelobes.  The  radar  PRF  is  changed  from  1000  pps  to  1150  pps  on  alternate  CPIs  to 
avoid  the  blind-speed  problem. 


FROM 
A/O 
10  8ITS 
2 6-MHz 
RATE 


760 

RANGE 

CELLS 

PER 

SWEEP 


Kig  14  — MTD  signal  proccw'r 


■R  J Me Aul.ty.  Tech.  Note  1972-14,  Lincoln  l aboratory.  Mass.  Inct  of  Tech  . 1972  , 

tR  M O'Donnell.  ( L Muche.  M Lahilt.  W II.  Drury,  and  L.  Cartledge.  h ASCON  Convention  Record  71-75.  1974 
*(  K.  Muche.  L Carlledge.  W.  It.  Drury.  H.  M.  Ilofstctter.  M l.ahitl.  P B McCorison.anJ  V.  J STerrino,  Proc.  ILEE 
A 2.  716.724.  (1974) 
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The  2.9  x 106  range-azimuth-doppler  cells  (760  x 360/0.75  x 8)  are  individually  thres- 
holded.  In  this  process  a clutter  map  is  generated  by  weighting  the  radar  return  in  the  zero- 
doppler  filter  over  the  last  eight  scans  (32  s)  using  a digital  filter.  Thus  tangential  targets  hav- 
ing zero  doppler  can  be  detected  if  the  target  level  exceeds  the  clutter-map  level  by  a specified 
constant.  That  is,  tangential  targets  can  be  detected  in  spotty  ground  clutter  by  using  the  prin- 
ciple of  in terclutter  visibility*.  The  thresholds  for  filters  2 through  6 are  set  using  a mean-level 
threshold.  Specifically  the  threshold  for  a given-number  filter  is  basedt  on  the  average  return 
in  the  given-number  filter  from  the  range  cells  ^1/2  n.mi.  (eight  cells)  on  either  side  of  the 
test  cell.  Since  clutter  spills  over  into  fillers  1 and  7,  two  thresholds  are  generated  for  these 
filters.  One  threshold  is  based  on  the  map,  a second  threshold  is  based  on  the  mean  level  over 
a range  interval,  and  the  higher  of  the  two  thresholds  is  used. 

The  MTD  represents  a great  improvement  in  signal  processing  for  FAA  air-surveillance 
radars.  A good  match  of  processor  to  radar  has  been  designed,  and  component  technology  has 
made  the  processing  practical  to  implement.  Presently,  a second-generation  MTD  is  being 
designed.  This  MTD  uses  no  MTI,  but  rather  each  filter  is  optimized  to  obtain  the  maximum 
signal-to-clutter-plus-noise  ratio  for  an  assumed  clutter  spectrum. 


Noncoherent  Moving-Target  Indicators 

Noncoherent  MTIs  are  described  in  Skolnik's  Introduction  to  Radar  Systems t and  Radar 
Handbook .§  They  differ  from  coherent  MTI  by  not  using  an  internal  coherent  reference  source 
but  rather  mixing  the  received  signal  with  itself.  Thus,  when  both  clutter  and  a target  are 
present,  the  beat  between  them  yields  a return  at  the  target  doppler.  On  the  other  hand,  when 
only  a target  is  present,  the  signal  return  is  at  zero  doppler  and  cannot  be  delected.  Conse- 
quently, for  noncoherent  MTI  to  be  useful,  gating  circuitry  is  required  for  passing  the  non- 
coherent MTI  output  when  clutter  is  present  and  passing  the  regular  video  when  clutter  is  not 
present.  Generally  fringe  areas  cause  major  problems  for  the  gating  circuitry,  making  perfor- 
mance unacceptable 

A different  kind  of  noncoherent  MTI  has  been  made  possible  by  high-power  microwave 
sources  ft  Lewis  and  Cantrell**  propose  transmitting  a short  pulse  and  subtracting  successive 
noncoherent  pulses  This  is  similar  to  an  area  MTI  discussed  in  Introduction  to  Radar  Systems. $ 
except  that  the  short  pulse  enables  the  subtraction  to  be  made  on  a pulse-lo-scan  pulse  rather  I , 

than  a scan-lo-scan  basis  Thus,  with  a I ns  pulse  and  a PRF  of  200  pps,  all  moving  targets 
above  60  knots  can  he  detected,  that  is,  there  are  no  blind  speeds. 


*D  K Barton  anil  W vv  Shrader.  I \SC ON  Cnni  Record  294-297,  I9h9  j 

♦ |)ei.nl>  jhout  v.inous  thresholding  lechmyues  can  be  found  in  ihe  section  on  noncoherent  processing 
\1  l Skolnik.  Ititnului  t"f r t>>  Rulur  St  \.vwv.  McCirjw  Hill.  New  \ork 

I Skolnik.  editor  Radar  Handbook.  McCr*w-H»ll.  New  York.  1970  | 

tt  \ . |.  CiranatMcin.  P Sprangle.  M Herndon  K K Parker,  and  S P Schlestngcr,  J Applied  Ph>sic>*  4f>.  3800-380'' 

‘ 197') 

"II  1.  I css  is  and  0 M Cantrell  Short  Pul',*  Noncoherent  Mil  patent  application.  Navy  Case  t>0$72.  M<L.  Nov. 

1975 
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NONCOHERENT  DETECTION 

The  earliest  noncoherent  signal  processing  was  performed  by  radar  operators  using  visual 
inputs  from  PPls  and  A-scopes.  Although  operators  can  perform  this  detection  task  accurately, 
operators  are  easily  saturated  and  become  quickly  fatigued.  To  remedy  this  situation  and  to 
provide  quick  reaction  times,  automatic  detection  and  tracking  (ADT)  systems  have  become 
quite  popular  during  the  1970s.  The  statistical  framework  necessary  for  the  development  of 
ADT  was  introduced  to  the  radar  community  in  the  1940s  by  Marcum*,  and  later  Swerlingt 
extended  the  work  to  fluctuating  targets.  They  investigated  many  of  the  statistical  problems  as- 
sociated with  the  noncoherent  detection  of  targets  in  Rayleigh  noise.  Their  mosi  important 
result  was  the  generation  of  curves  of  probability  of  detection  ( PD ) versus  signal-to-noise  ratio 
(S/N)  for  a detector  which  sums  N enveloped  detected  samples  (either  linear  or  square  law) 
under  the  assumption  of  equal  signal  amplitudes.  However,  in  a search  radar,  as  the  beam 
sweeps  over  the  target,  the  returned  signal  amplitude  is  modulated  by  the  antenna  pattern. 
Many  authors  investigated  various  detectors  (weightings),  comparing  detection  performance 
and  angular  estimation  results  to  the  optimal  values.  The  detectors  investigated  included  the 
moving  window,  feedback  integrator,  two-pole  filter,  binary  integrator,  and  batch  processor. 

In  the  original  work  on  these  detectors,  the  environment  was  assumed  known  and  homo- 
geneous, so  that  fixed  thresholds  could  be  used.  However  a realistic  environment,  containing 
land,  sea,  and  rain  for  example,  will  cause  an  exorbitant  number  of  false  alarms  for  a fixed 
threshold  system.  Two  approaches,  adaptive  thresholding  and  nonparametric  detectors,  have 
been  used  to  solve  the  false-alarm  problem.  Both  solutions  are  based  on  the  assumption  that 
homogeneity  exists  in  a small  region  about  the  range  cell  that  is  being  tested.  The  adaptive 
thresholding  method  assumes  that  the  noise  density  is  known  except  for  a few  unknown 
parameters.  The  surrounding  reference  cells  are  then  used  to  estimate  the  unknown  parame- 
ters, and  a threshold  based  on  the  estimated  density  is  obtained.  Nonparametric  detectors  ob- 
tain a constant  false-alarm  rate  (CFAR)  by  ranking  the  test  sample  with  the  reference  cells. 
Under  the  hypothesis  that  all  the  samples  (test  and  reference)  are  independent  samples  from 
an  unknown  density  function,  the  test  sample  has  a uniform  density  function;  consequently  a 
threshold  which  yields  CFAR  can  be  set. 

Classical  Theory 

The  radar  detection  problem  is  a binary-hypothesis-testing  problem; 

//0:  no  target  present 
or 

H | : target  present. 

Many  criteria  can  be  used  to  solve  this  problem,  but  the  most  appropriate  for  radar  is  the 
Ncyman-Pearsont  criterion.  This  criterion  maximizes  Pp  for  a given  probability  of  false  alarm 


*J  I.  Marcum.  IRE  Tunv  Information  Theory  6.  S9-767  (196(1) 
tp  Sterling.  IRK  Trans.  Information  Theory  6,  269*308  (I960), 
tj  Ncyman  and  F.  S Pearson.  Biomctrika  20. \,  175-240,  263-294  M928) 
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(Pla)  by  comparing  the  likelihood  ratio  ( L ) to  an  appropriate  threshold  T. 
present  if 


L ) 


/>(*! -VJ  H\  ) 

p(x I x„\H0  ) 


A target  is  declared 


(45) 


where  p Cv, ,v„|  //,  ) and  p(xx . ....  ar„|  H0  ) are  the  joint  densities  of  the  n samples  under  the 

conditions  of  target  presence  and  target  absence  respectively.  For  a linear  envelope  detector 
and  white  Gaussian  noise  the  samples  have  a Rayleigh  density  under  H0  and  a Ricean  density 
under  //] , and  the  likelihood  detector  reduces  to 


n 


ru 


T. 


(46) 


where  /q  is  the  Bessel  function  of  zero  order.  For  equal-amplitude  (A.  — .4)  small  signal 
pulses  (A,  < < o-),  the  detector  reduces  to  the  square-law  detector: 

t x}  > T.  (47) 

i -i 


This  detector  and  the  linear  detector  were  first  studied  by  Marcum*  and  were  studied  in 
succeeding  years  by  numerous  people.  The  most  important  facts  concerning  these  detectors 
are  the  following: 


• The  detection  performances  of  the  linear  and  square-law  detectors  are  similar  and  are 
close  to  the  performance  of  the  optimal  detector.* 


• Since  the  signal  return  of  a scanning  radar  is  modulated  by  the  antenna  pattern,  only 
0.84  of  the  pulses  between  the  half-power  points  should  be  integrated,  and  the  antenna  beam- 
shape  factor  (ABSF)  is  1.6  dB.f  The  ABSF  is  the  number  by  which  the  midbeam  S/N  must  be 
reduced  so  that  the  detection  curves  generated  for  equal  signal  amplitudes  can  be  used  for  the 
scanning  radar. 

• The  collapsing  loss  for  the  linear  integrator  can  be  much  greater  than  the  loss  for  a 
square-law  integrator.*  The  collapsing  loss  is  the  additional  signal  required  to  maintain  the  same 
PD  and  P/a  when  unwanted  noise  samples  along  with  the  desired  signal-plus-noise  samples  are 
integrated. 

Most  signal  processors  are  required  not  only  to  detect  targets  but  to  make  angular  esti- 
mates of  the  azimuth  position  of  the  target.  Swerling§  calculated  the  standard  deviation  ot  the 
optimal  estimate  by  using  the  Cramer-Rao  lower  bound.  The  results  are  shown  in  Fig.  15, 
where  a normalized  standard  deviation  is  plotted  against  S/N  per  pulse.  This  result  holds  for  a 
moderate  or  large  number  of  pulses  integrated,  and  the  optimal  estimate  involves  finding  the 
location  where  the  correlation  of  the  returned  signal  and  the  derivative  of  the  antenna  pattern 
is  zero.  Although  this  estimate  is  rarely  implemented,  its  performance  is  approached  by  simple 


i l 


"J  I Marcum.  I K I Trans  Information  Theory  6.  s').jfi7  (I960). 
+L  V.  Blake.  Proc.  IRE  41. 770-774  (1953). 

>Ci  V.  Trunk,  Proc  IEEE  60.  743-744  (1973). 
ljl>.  Sterling.  Proc  IRE  44.  1146-1155  (1956). 
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Fig.  15  — Comparison  of  angular  esli- 
maies  wiih  ihe  Cramcr-Rao  lower 
bound.  In  ihe  ordinate  expression,  <r  is 
the  standard  deviation  of  Ihe  estimation 
error  and  N is  the  number  of  pulses 
within  the  3-dB  bcamwidth.  which  is 
20- 


estimates,  such  as  the  maximum-value  and  threshold-crossing  procedures,  as  can  be  seen  in 
Fig.  15. 


Integrators 


Almost  all  signal  processors  use  linear  rather  than  square-law  detectors,  since  a linear 
detector  is  easily  built  by  using  a matched  filter  and  a half-wave  rectifier  followed  by  a low-pass 
filter.  However  many  different  integrators  are  used  to  accumulate  the  linear-envelope-detected 
pulses.  A few  of  the  most  common  integrators  are  shown  in  Fig.  16.  Some  advantages  and 
disadvantages  of  these  integrators  are  as  follows. *ft 


Moving  window 

The  moving  window  performs  a running  sum  of  /V pulses;  as  the  latest  pulse  is  added  to 
the  sum,  the  pulse  that  is  V PR Fs  in  the  past  is  subtracted  from  the  sum  The  detection  per- 
formance of  this  detector  is  only  0.5  dB  worse  than  the  optimal  detector  which  weights  the  re- 
turned signal  by  the  fourth  power  of  the  voltage  antenna  pattern.  The  angular  estimate  is  ob- 


M)  S P.ilmer  and  I)  C\  Cooper.  ILL l*.  Tran%  Information  Theory  IT-10.  2%-302  U%4). 
tCi  \1  Dillard.  ILLL  Trans  Information  Theory  IT-13,  2-6  (l%7). 

*B  H Cantrell  and  G.  V.  Trunk.  1 1 LI  Trans  Aerospace  and  Ltectronic  Systems  \ES-0.  f>49-t»53  U97a) 
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Fig.  16  ““  Common  integrators 


tained  by  either  taking  the  maximum  value  of  the  running  sum  or  taking  the  midpoint 
between  the  first  and  last  crossing  of  the  detection  threshold.  Both  methods  have  a bias  of  N/2 
pulses  which  is  easily  corrected.  The  standard  deviation  of  the  estimation  error  of  both  estima- 
tors is  about  20%  higher  than  the  Cramer-Rao  lower  bound.  The  major  disadvantage  of  this 
detector  is  that  the  last  N pulses  for  each  range  cell  must  be  saved.  For  radars  with  large 
beamwidths  and  thus  many  pulses,  the  moving  window  requires  extensive  hardware.  However 
with  the  lower  cost  and  size  of  memory  this  disadvantage  is  rapidly  disappearing. 


Feedback  integrator 

The  amount  of  storage  required  can  be  reduced  significantly  by  using  a feedback  integra- 
tor. which  requires  the  storage  of  only  one  number.  Although  the  feedback  integrator  applies 
an  exponential  weighting  into  the  past,  its  detection  performance  is  only  I dB  less  than  the  op- 
timal integrator.  Unfortunately  difficulties  are  encountered  when  using  the  feedback  integrator 
to  estimate  the  azimuth  position  The  threshold-crossing  procedure  yields  estimates  only  20% 
create  than  the  lower  bound,  but  the  bias  is  a function  of  S/N  and  must  be  estimated.  On  the 
other  hand  the  maximum  value,  although  having  a constant  bias,  has  estimates  which  are  100% 
greater  than  the  lower  bound  This  author's  opinion  is  that  this  detector  has  limited  utility. 
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T wo-pole  filter 

The  two-pole  filter  requires  the  storage  of  an  intermediate  calculation  in  addition  to  the 
integrated  output.  However  with  this  rather  simple  device  a weighting  pattern  similar  to  the 
antenna  pattern  can  be  obtained;  consequently  good  performance  would  be  expected.  The 
detection  performance  is  within  0.15  dB  of  the  optimal  detector,  and  its  angular  estimates  are 
about  20%  greater  than  the  Cramer-Rao  lower  bound.  If  the  desired  number  of  pulses  in- 
tegrated is  changed  (because  of  change  in  rotation  of  the  radar  or  use  of  another  radar),  it  is 
necessary  to  change  only  the  feedback  values  A',  and  K2.  Their  optimal  values  are  set  by 

„ , — linijr/v  I —ih2  , , 

A|  = 2 r cos  ( oi (/T ) 

and 

,,  _ -2.'/c.vr/x/l  — 2 

A2  " £ 

where  <b  = 0.63,  = 2.2,  and  N is  the  number  of  pulses  between  the  3-dB  points  of  the 

antenna. 

Binary  Integrator 

The  binary  integrator  is  also  known  as  the  dual-threshold  detector,  M-out-of-A  detector, 
or  rank  detector.  The  input  samples  are  quantized  to  0 or  1 depending  on  whether  or  not  they 
are  less  than  a threshold  7"| . The  last  N zeros  and  ones  are  summed  and  compared  to  a second 
(detection)  threshold  T2  = M.  The  detection  performance  of  this  detector  is  2 dB  less  than 
the  moving-window  integrator  because  of  the  hard  limiting  of  the  data,  and  the  angular  estima- 
tion error  is  25%  greater  than  the  Cramer-Rao  lower  bound.  This  detector  is  used  because  it  is 
easily  implemented,  it  ignores  interference  spikes  which  cause  trouble  with  integrators  that 
directly  use  signal  amplitude,  and  it  works  extremely  well  when*t  the  noise  has  a non- 
Rayleigh  density. 

A comparison  of  the  binary  integrator  (three  out  of  three),  the  median  detector  (two  out 
of  three),  and  the  mean  detector  (moving  window)  in  log-normal  interference  is  shown  in  Fig. 
17.  The  optimal  binary  integrator  is  much  better  than  straightforward  integration.  The  optimal 
values  for  the  second  threshold  were  found  by  Schwartzi  for  Rayleigh  interference  and  by 
Schleher§  for  log-normal  interference. 

Batch  Processor 

The  batch  processor  is  used  when  there  are  a large  number  of  pulses  in  the  3-dB 
bcamwidth.  If  AW  pulses  are  in  the  3-dB  beamwidth.  A'  pulses  are  summed  and  either  a 0 or  1 
is  declared  depending  on  whether  or  not  the  sum  is  less  than  a threshold  T{.  The  last  A' zeros 
and  ones  are  summed  and  compared  to  a second  threshold  \f. 


(48) 

(49) 


*1)  C Schleher.  Il  f I 1975  International  Kad.tr  Conf . 2b2-2b7.  1975 

Hi  V.  Trunk.  Non-Raclcigh  Sea  Clutter  I'ropcrt ios  ,nul  Detection  of  Targets.  NRl.  Report  79Kb,  June  |97b. 
»M  Sv.hrc.irt/  II  I L Traits  Information  T heor>  2.  135-139  (195b) 

71)  ('  Sehleh  -r  II  I I I97'  International  R.nJar  Conf.  2b2-2h7.  1975 
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Fig.  17  — ComparKon  of  various  ilcleclors  in  log-normal 
imcrlcrcncc  (,V  = 3, /%  = I0'6) 


I The  batch  processor,  like  the  binary  integrator,  is  easily  implemented,  ignores  interfer- 

ence spikes,  and  works  extremely  well  when  the  noise  has  a non-Rayleigh  density,  but  further- 
more in  comparison  with  the  binary  integrator  the  batch  processor  requires  less  storage,  detects 
better  (less  than  2 dB  from  moving  window),  and  estimates  angles  more  accurately. 


The  batch  proccessor  has  been  implemented  by  the  Applied  Physics  Laboratory*  of 
Johns  Hopkins  University  with  great  success.  To  obtain  a more  accurate  azimuth  estimate, 
they  use 


a 


£4 4 

14 


(50) 


where  A/  are  the  amplitudes  of  the  sums  greater  than  Tt  and  0,  are  the  corresponding  antenna 
azimuth  angles.  When  many  pulses  are  on  target  ( N > 20),  this  detector  is  generally  favored 
by  this  author. 


False  Alarms 

If  fixed  thresholds  are  used  with  the  previously  discussed  integrators,  the  detectors  will 
saturate  the  tracking  computer  associated  with  the  system  and  disrupt  the  system.  Three  im- 
portant facts  should  be  remembered: 

• It  makes  little  sense  to  have  an  automatic  detection  system  without  an  associated 

tracking  system; 

• The  sensitivity  of  the  detector  should  be  as  high  as  possible  without  saturating  the 

tracking  computer; 


l 


I 
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• False  alarms  and  false  targets  are  not  a problem  if  they  are  removed  by  the  tracking 
computer.  Tracking  (scan-to-scan  processing)  is  the  only  way  to  remove  stationary 
point  clutter  or  target  MTI  residues. 

One  can  reduce  the  number  of  false  alarms  with  a fixed-threshold  system  by  setting  a 
high  threshold,  but  this  would  reduce  sensitivity  in  regions  of  low-noise  (clutter)  return.  A 
detector  is  required  which  will  detect  a target  when  it  has  a higher  return  than  its  immediate 
background.  Two  such  types  of  detectors  are  adaptive-thresholding  and  nonparametric  detec- 
tors. Both  of  these  detectors  assume  that  the  samples  in  the  range  cells  surrounding  the  test 
cell  (called  reference  or  neighboring  cells)  are  independent  and  identically  distributed;  further- 
more it  is  usually  assumed  that  the  time  samples  are  independent.  Both  detectors  test  whether 
the  test  cell  has  a return  sufficiently  larger  than  the  reference  cells.  A survey  of  CFAR  pro- 
cedures can  be  found  in  Hansen*. 

Adaptive  Thresholding 

The  basic  assumption  of  the  adaptive-thresholding  technique  is  that  the  noise  density  is 
known  except  for  a few  unknown  parameters.  The  surrounding  reference  cells  are  used  to  es- 
timate the  unknown  parameters,  and  a threshold  based  on  the  estimated  density  is  then  ob- 
tained. The  simplest  adaptive  detector  is  the  cell-averaging  CFAR  investigated  by  Finn  and 
Johnsont.  If  the  noise  has  a Rayleigh  density,  only  the  parameter  ir  needs  to  be  estimated, 
since  the  mean  of  a Rayleigh  distribution  is  irjrrji  and  the  variance  is  <r  2(2v— tt/2  ).Thus,  by 
estimating  the  mean,  one  obtains  an  estimate  <r  which  can  be  used  to  set  a threshold  T to  yield 
the  desired  Pla.  However,  since  T is  set  by  an  estimate  <r,  it  must  be  slightly  larger  than  the 
threshold  one  would  use  if  <r  were  known  a priori.  The  raised  threshold  causes  a loss  in  target 
sensitivity  and  is  referred  to  as  a CFAR  loss.  This  loss  has  been  calculated  by  Mitchell  and 
Walker*,  and  some  results  are  summarized  in  Table  2.  As  can  be  seen,  for  a small  number  of 
reference  cells,  the  loss  is  large  because  of  the  poor  estimate  of  <r. 

This  thresholding  technique  is  more  effective  in  maintaining  CFAR  when  it  is  applied  to 
the  binary  integrator  or  batch  processor,  as  shown  in  Fig.  18.  This  is  because  when  the 
number  of  pulses  integrated  by  the  binary  integrator  is  moderate,  the  P^t,  on  a single  pulse  is 
rather  large;  for  example  P/a  =0.1  for  a single  pulse  yields  Pfo  = 10  5 for  a seven-out-of- 
ten  integrator.  Thus,  since  most  non-Rayleigh  densities  are  Rayleigh-like  to  the  10th  percen- 
tile, this  type  of  processor  will  maintain  a low  Pfo  in  most  non-Rayleigh  environments.  This 
demonstrates  a general  rule,  to  maintain  a low  Pfa  in  various  environments,  adaptive  thres- 
holding should  be  placed  in  front  of  the  integrator.  For  any  noise  distribution.  CFAR  can  be 
maintained  by  counting  the  number  of  ones  out  of  the  comparator  per  scan  and  using  this 
number  to  control  K:  that  is,  if  the  number  is  too  large,  K is  increased. 

Front-end  thresholding,  which  maintains  amplitude  information  by  dividing  the  average 
reference  value  into  the  test  cell,  was  investigated  by  Hansen  and  \Vard§  and  is  shown  in  Fig. 
19.  This  type  of  processing  is  especially  effective  when  there  is  strong  interference  which  is 
variable  on  a pulse-io-pulse  basis. 
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When  the  noise  has  a non-Rayleigh  density,  such  as  the  chi-square  density  or  log-normal 
density,  two  parameters  must  be  estimated,  and  the  adaptive  detector  is  more  complicated  If 
several  pulses  are  integrated  with  any  of  the  amplitude  integrators,  the  integrated  output  will 
be  approximately  Gaussian  distributed.  Then  the  two  parameters  which  must  be  estimated  are 
the  mean  and  the  variance.  These  estimates  are  given  by 

(51) 

I 

and 

r?  2 = L - *2  (52) 

where  the  summation  is  over  the  /V  range  cells  surrounding  the  test  cell. 

When  successive  pulses  in  the  same  range  cell  are  correlated  (as  with  returns  from  rain 
or  sea  clutter),  many  false  alarms  will  occur  if  only  the  mean  value  (51)  is  estimated.  A thres- 
hold of  the  form 

T = X + K<r  (53) 

will  provide  a low  P/a  for  the  amplitude  integrators:  moving  window,  feedback  integrator,  and 
two-pole  filter.  Nothing  can  be  done  to  the  binary  integrator  to  yield  a low  Pta  in  correlated 
noise;  thus  it  should  not  be  used  in  this  situation.  On  the  other  hand,  if  the  correlation  time  is 
less  than  a batching  interval,  the  batch  processor  will  yield  a low  P/a  without  modifications. 

Nonparametnc  Detectors 

The  most  common  way  nonparametric  detectors  obtain  CFAR  is  by  ranking  the  test  sam- 
ple with  the  reference  cells.  Under  the  hypothesis  that  all  the  samples  are  independent  sam- 
ples from  an  unknown  density  function,  the  test  sample  has  a uniform  density  function.  For 
instance,  with  reference  to  the  rank  detector  in  Fig.  20,  the  test  cell  is  compared  to  15  of  its 
neighbors.  Since  in  the  set  of  16  samples  the  test  sample  has  equal  probability  of  being  the 
smallest  sample  (rank  = 0 or  equivalently  any  other  rank),  the  probability  that  the  test  sample 

takes  on  values  0,  1 15  is  1/16.  A simple  rank  detector*  can  be  constructed  by  comparing 

the  rank  (number  of  reference  cells  that  the  test  cell  exceeds)  to  a threshold  K,  and  the  output 
is  1 if  the  rank  is  larger  and  0 otherwise.  The  zeros  and  ones  are  summed  in  a moving  win- 
dow. This  detector  incurs  a CFAR  loss  of  about  2 dB  and  is  extremely  effective,  if  the  time 
samples  are  independent.  Only  certain  values  of  P/a  can  be  obtained.  Thus,  if  the  number  of 
pulses  integrated  is  small,  low  P.a  values  cannot  be  obtained. 

If  the  time  samples  are  dependent,  the  rank  detector  will  not  yield  CFAR.  A modified 
rank  detector,  called  the  modified  generalized  sign  testf  (MGST)  is  an  attempt  to  maintain  a 
low  P/a  and  is  that  shown  in  Fig.  20.  This  detector  can  be  divided  into  three  parts:  a ranker, 
an  integrator  (in  this  case  a two-pole  filter),  and  a thresholding  device.  A target  is  declared 
when  the  integrated  output  exceeds  two  thresholds.  The  first  threshold  is  fixed  (equals 
M + T|/A'  from  Fig.  20)  and  yields  CFAR  when  the  reference  cells  are  independent  and 
identically  distributed.  The  second  threshold  is  adaptive  and  maintains  a low  P,a  when  the 

*V  (i  Manser)  and  B A Olsen.  IEEE.  Trans  Aerospace  and  Electronic  Ssstems  4.  942-950  (1971) 
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reference  samples  are  correlated.  The  device  uses  the  mean-deviate  estimate,  where  ex- 
traneous targets  in  the  reference  cells  have  been  excluded  from  the  estimate  by  use  of  a pre- 
liminary threshold  T2 , to  estimate  the  standard  deviation  of  the  correlated  samples. 

The  rank  and  MGST  detectors  are  basically  two-sample  detectors.  They  decide  a target  is 
present  if  the  ranks  of  the  test  cell  are  significantly  greater  than  the  ranks  of  the  reference 
cells.  Target  suppression  occurs  at  all  interfaces  (such  as  a land-and-sea  interface),  where  the 
homogeneity  assumption  is  violated.  However,  some  tests  exist  (Hansen*  investigated  the 
Spearman  Rho  and  Kendall  Tau  tests)  which  depend  on  only  the  test  cell.  These  tests  work  on 
the  fact  that,  as  the  antenna  beam  sweeps  over  a point  target,  the  signal  return  increases  and 
then  decreases.  Thus  for  the  test  cell  the  ranks  should  follow  a pattern  first  increasing  and 
then  decreasing.  Although  these  detectors  do  not  require  reference  cells  and  hence  have  the 
useful  property  of  not  requiring  homogeneity,  these  detectors  are  not  generally  used  because  of 
the  large  CFAR  loss  taken  for  moderate  sample  sizes:  for  N - 16  the  loss  is  10  dB,  and  for 
N = 32  the  loss  is  6 dB. 

The  paper  by  Hansen*  is  worth  noting  because  it  introduced  the  concept  of  importance- 
sampling for  calculation  of  false-alarm  thresholds.  The  fundamental  principle  of  the 
importance-sampling  technique  is  to  modify  the  probabilities  that  govern  the  outcome  of  the 
basic  experiment  of  the  simulation  in  such  a way  that  the  event  of  interest  (the  false  alarm) 
occurs  more  frequently.  This  distortion  is  then  compensated  for  by  weighting  each  event  by 
the  ratio  of  the  probability  that  this  specific  event  would  have  occurred  if  the  true  probabilities 
had  been  used  in  the  simulation  to  the  probability  that  this  event  would  occur  with  the  distort- 
ed probabilities.  Consequently  by  proper  choice  of  the  distorted  probabilities  the  number  of 
repetitions  can  be  reduced  greatly.  Further  details  on  importance  sampling  can  be  found  in 
Trunk  et  al.f,  Hansent,  and  Hillier  and  Licberman§. 

In  summary,  when  only  a small  number  of  returns  are  available  (less  than  eight),  ampli- 
tude information  must  be  used,  and  this  author  favors  the  moving-window  integrator.  When  a 
moderate  number  (between  eight  and  20)  are  available,  a rank  detector  should  be  used  if  sam- 
ples are  independent;  and  a two-pole  filter  with  thresholding  of  the  form  T = X 4-  Kir  should 
be  used  if  the  samples  are  dependent.  If  a larae  number  of  pulses  (greater  than  20)  are  avail- 
able, the  batch  processor  or  MGST  processor  should  be  used.  These  rules  should  serve  only  as 
a general  guideline.  It  is  highly  recommended  that  a sample  of  the  radar  environment  be  col- 
lected and  analyzed  and  that  various  debtors  be  simulated  on  a computer  and  tested  against 
recorded  data. 

Sequential  Detectors 

Sequential  detectors,  which  can  be  used  with  phased-array  radars,  are  based  on  the  idea 
that  in  many  cases,  depending  on  the  returned  samples,  a decision  can  be  made  on  a few  sant- 
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plex.  The  sequential  likelihood-ratio  test  (SLRT)  works  as  follows:  given  independent  samples 
A,  v„r  calculate  the  likelihood  ratio 


L 


in 


P<x ! *W|  H\  > 

Pi*  1 X»,\H0) 


(54) 


It  Lm  > A,  accept  (target  present);  if  Lnl  ^ B.  accept  H0  (no  target  present);  and  if 
® < A,  take  another  sample.  The  SLRT  has  the  useful  properties  that  the  thresholds 

are  set  by  the  simple  formulas  A = PD/P,a  and  B = (1  - Pn  )/(l  - P,u)  and  that  for  all 
tests  with  a given  PD  and  Pfa  the  SLRT  requires  the  smallest  average  sample  size.  Further  de- 
tails about  the  SLRT  can  be  found  in  Lindgren*. 


An  early  application  of  sequential  detection  to  radar  was  discussed  by  Marcus  and  Swer- 
lingf.  Unfortunately,  in  radar  the  application  of  SLRT  is  not  straightforward,  since  one  is  re- 
quired to  make  a decision  in  every  range  cell  before  the  test  can  be  ended  and  the  agile  beam 
moved.  The  modified  problem  considered  was 


//0:  noise  present  in  all  range  cells 


or 


//| : exactly  one  signal  present  in  the  rth  range  cell  (i  unknown). 

They  performed  some  numerical  calculations  and  concluded  that: 

• The  greatest  savings  in  average  sample  size  comes  when  no  signal  is  present  (//„ 

true); 

• In  comparison  with  a fixed-sample-size  test,  SLRT  provides  a greater  savings  when 

the  number  of  range  cells  is  small  and  when  S/N  is  small; 

• It  is  not  necessary  to  truncate  the  test. 

TRACKING  SYSTEM 

Track-while-scan  systems  (tracking  systems  for  surveillance  radars  whose  nominal  scan 
time  if  from  4 to  12  s)  will  now  be  considered.  If  the  probability  of  detection  ( Pt) ) per  scan  is 
high,  if  accurate  measurements  are  made,  if  the  target  density  is  low,  and  if  there  are  few  false 
detections  (crossings  of  a threshold,  with  no  judgment  being  made  on  whether  or  not  it  be- 
longs to  a valid  target),  the  design  of  the  correlation  logic  and  tracking  filter  is  straightforward. 
However  in  a realistic  radar  environment  these  assumptions  are  never  valid,  and  the  design 
problem  is  complicated. 

White  and  Silbermant  list  many  problems  encountered  in  actual  situations.  Among  these 
problems  are  target  fades  (due  to  multipath  propagation,  clutter  masking,  interference,  blind 


*H  W t indgren.  Sr(iii\iuul  Theory,  MacMillan,  New  York,  l%2 
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speeds,  and  atmospheric  conditions),  false  alarms  (due  to  noise,  clutter,  interference,  and  jam- 
ming), and  poor  radar  parameter  estimates  (due  to  noise,  unstabilized  radar  platforms,  un- 
resolved targets,  target  splits  (two  detections  for  a single  target),  multipath  propagation,  and 
propagation  effects). 

A general  outline  of  a track-while-scan  system  will  be  considered  first.  Then  the  tracking 
filer,  maneuver-following  logic,  track  initiation,  and  correlation  logic  will  be  discussed  in  detail. 
Finally,  methods  of  integrating  data  from  several  radars  will  be  discussed. 

System  Outline 

Almost  all  track-while-scan  systems  operate  on  a sector  basis.  A typical  series  of  opera- 
tions is  shown  in  Fig.  21.  For  instance,  if  the  radar  has  reported  all  the  detections  in  sector  1 1 
and  is  now  in  sector  12,  the  tracking  program  would  start  by  correlating  (trying  to  associate) 
the  clutter  points  (stationary  tracks)  in  sector  10  with  detections  in  sectors  9,  10,  and  11. 
Those  detections  that  are  associated  with  clutter  points  are  deleted  (are  not  used  for  further 
correlations)  from  the  detection  file  and  are  used  to  update  the  clutter  points.  Updating  clutter 
points  usually  implies  replacing  the  old  point  by  the  associated  detection. 


Fig  21  — Various  opcraiions  of  a 

Imck-whilc-scan  system  performed  on 
a sector  basis 


Next,  firm  tracks  in  sector  8 are  correlated  with  detections  in  sectors  7,  8,  and  9.  By  this 
time  all  clutter  points  have  been  removed  from  sectors  9 and  below.  Those  detections  which 
are  associated  with  firm  tracks  are  deleted  from  the  detection  file  and  are  used  to  update  the 
appropriate  track.  The  filter  for  performing  this  updating  will  be  described  in  the  next  subsec- 
tion. 


Usually,  some  provision  is  made  for  giving  perference  to  firm  tracks  (with  respect  to  ten- 
tative tracks)  in  the  correlation  process.  By  performing  the  correlation  process  two  sectors 
behind  firm  track  correlations  (Fig.  21),  it  is  impossible  for  tentative  tracks  to  steal  detections 
belonging  to  firm  tracks.  In  other  tracking  systems  the  correlation  for  firm  and  tentative  tracks 
is  performed  in  the  same  sector;  however  the  generalized  distance  D between  tracks  and  detec- 
tions is  incremented  by  ID  if  the  track  is  tentative. 

Finally  detections  which  are  not  associated  with  either  clutter  points  or  tracks  are  used 
for  initiation  purposes.  The  most  common  initiation  procedure  is  to  initiate  a tentative  track; 
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I, iter  the  track  is  dropped  or  else  made  a firm  track  or  clutter  point.  Cantrell  et  ul.*  suggest  that 
both  a clutter  point  and  tentative  track  be  established.  If  the  detection  came  from  a stationary 
target,  the  clutter  point  will  be  updated  and  the  tentative  track  will  eventually  be  dropped.  On 
the  other  hand,  if  the  detection  came  from  a moving  target,  the  tentative  tracks  will  be  made 
firm  and  the  clutter  point  will  be  dropped.  The  latter  method  requires  less  computer  computa- 
tion time  when  most  of  the  detections  are  clutter  residues. 

The  correlation  procedure  is  made  in  a sector  framework  to  avoid  the  necessity  of  corre- 
lating all  tracks  with  all  detections.  The  procedure  can  be  implemented  very  easily  by  defining 
two  computer  arrays:  a sector  file  and  a track  file.  The  sector  file  for  sector  I contains  the  first 
track  number  in  sector  I,  and  the  track  file  for  track  J contains  the  next  track  number  in  the 
same  sector  as  track  J or  a zero,  indicating  that  the  track  is  the  last  track  in  the  sector.  Further 
information  about  the  details  of  a tracking  system  can  be  found  in  Cantrell  et  al.t,  Wilson  and 
Cantrell*,  and  Trunk  et  al§. 

Tracking  Filters 

Before  proceeding,  the  coordinate  system  in  which  the  tracking  will  be  performed  will  be 
discussed.  The  quantities  measured  by  the  radar  are  spherical:  range,  azimuth,  elevation,  and 
possibly  range  rate.  Thus  it  may  seem  natural  to  perform  tracking  in  spherical  coordinates. 
However  this  causes  difficulties,  since  motion  of  constant-velocity  targets  (straight  lines)  will 
cause  acceleration  terms  in  all  coordinates.  A simple  solution  to  this  problem  is  to  track  in  a 
Cartesian  coordinate  system.  While  it  may  appear  that  the  appropriate  transformations 
(,v  = R cos  Ut,  cos  H ir  etc.  where  R is  the  range,  fl(,  is  the  elevation  angle,  and  Ha  is  the  az- 
imuth angle)  will  destroy  the  accurate  range  track,  Cantrell#  has  shown  that  the  inherent  ac- 
curacy is  maintained.  Quigiey  and  Holmes**  note  that  maneuvering  targets  cause  a large  range 
error  but  a rather  insignificant  azimuth  error  and  thus  suggest  using  a target-oriented  Cartesian 
coordinate  system.  Specifically,  the  x axis  is  taken  along  the  azimuth  direction  of  the  target 
and  the  r axis  is  taken  in  the  cross  range  direction. 

Skalanskytt  performed  one  of  the  first  analyses  of  the  tracking  filter  for  a track-while- 
scan  system.  He  considered  the  n —fi  filter  described  by 

*,(*)  = x p (A)  + «[*„,  (A)  - xn  (A ) ],  (55) 

(A)  = Vs  (k  - 1 ) + H lx*  (*)  - xp  (k)]/T,  (56) 
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A p(k  + 1)  = *v  (A)  + vs  (k)T.  (57) 

where  xs  (A)  is  the  smoothed  position,  Vs  (A)  is  the  smoothed  velocity,  xp  (A)  is  the  predicted 
position,  a„,  (A)  is  the  measured  position,  T is  the  scanning  period  (time  between  detections), 
and  <*  and  fi  are  the  system  gains. 


The  optimal  filter  for  performing  the  tracking  when  the  equation  of  motion  is  known  is 
the  Kalman  filter,  first  discussed  by  Kalman*  and  later  by  Kalman  and  Bucyf.  The  Kalman 
filter  is  a recursive  filter  which  minimizes  the  least-square  error.  The  state  equation  in  xy  coor- 
dinates for  a constant-velocity  target  ist 


A'O  + l ) = <b(t)  XU)  + 10)  AU), 


where 


aO) 

1 

T 

0 

0 

- 

A'O)  = 

x(t) 

yU) 

• *<'>  -S 
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0 

0 

1 

0 

T 

>'0)l 

0 

0 

0 

1 

I'  0)  = 


and  A 0)  — 


av  0)  ‘ 


with  A'O)  being  the  state  vector  at  time  / consisting  of  position  and  velocity  components  a 0), 
a'O),  yU),  and  TO);  / + 1 being  the  next  observation  time;  T being  the  lime  between  obser- 
vations; and  ax  0)  and  av(t)  being  random  accelerations  with  covariance  matrix  Q(r).  The 
observation  equation  is 


where 


TO)  = 


TO)  = A/0)  XU)  + VU). 


£<;>■"<'»  -ii  s ? si-"*  *,(')  - 


Vv0)  • 


with  TO)  being  the  measurement  at  time  / consisting  of  positions  xm  (/)  and  >'„,(/)  and  ^0) 
being  a zero-mean  noise  whose  covariance  matrix  is  RU). 


The  problem  is  solved  recursively  by  first  assuming  the  problem  is  solved  at  time  / —1. 
Specifically,  it  is  assumed  that  the  best  estimate  X(i  — 1 [ z — 1 ) at  time  / —1  and  its  error  covari- 
ance matrix  PU  — 1 1 Z — 1 ) are  known,  where  the  caret  in  the  expression  of  the  form  A'0|s) 
signifies  an  estimate  and  the  overall  expression  signifies  A'O)  is  being  estimated  with  observa- 
tions up  to  T(s).  The  six  steps  involved  in  the  recursive  algorithm  are 

1.  Calculate  the  one-step  prediction 

V0|/-1)  =».U/-1  )A'0— 1 1 /-l);  (60) 


’H  I Kjlm.in.  t r.ms  ASM K Scries  0.  I Basic  Knginoering  82.  55-45  (I960). 

fit  t Kalman  and  K S Buc> , Trans  \S\tl.  Series  It.  J Basic  I ngincering  St  95-10*  (1961). 

H R Castclla  and  t;.  Ci  Dunnehacke.  Ill  I Trans  Aerospace  and  I leclrnnic  Syslems  10,  89|.g9s  (1974) 
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2.  Calculate  the  covariance  matrix  for  the  one-step  prediction 

P(r\  t -1  ) = <5  (/  -1  )PU  -1 1 r -1  )<b  T(r  -1  ) + I (/  -1  )Q(r  -1 )!'  T(t  -1 );  (61) 

3.  Calculate  the  predicted  observation 

Y (r| t -1 ) = M(i)X(i\t  -1 ):  (62) 

4.  Calculate  the  filter  gain 

A (/)  = PU\t  -1  ) MrU)  [MU)PU\t  -1  )MtU)  + R(r)]~':  (63) 

5.  Calculate  the  new  smoothed  estimate 

XU\t)  = XU\t-\)  + A (/)  [>(r)  - YU\r  -1)J;  (64) 

6.  Calculate  the  new  covariance  matrix 

PU\t)  = \I  - SU)MU)]P(r\r-\).  (65) 


In  Summary,  with  an  estimate  XU—  l|r—  1)  and  its  covariance  matrix  PU—\\t—  1)  as  the 
start,  after  a new  observation  YU)  is  received  and  the  six  quantities  in  the  recursive  algorithm 
are  calculated,  a new  estimate  X(r\r)  and  its  covariance  matrix  PU\t)  are  obtained. 


It  is  fairly  simple  to  show  that  for  a zero  random  acceleration.  Q(i)  = 0 and  a constant 
measurement  covariance  matrix  R(i)  = R , the  o —j3  filter  can  be  made  equivalent  to  the  Kal- 
man filter  by  setting 


2 (2  k - 1 ) 
k (k  + 1 ) 


(66) 


and 


on  the  Ath  scan* 


P 


6 

k (A  + 1 ) 


(67) 


Thus  as  time  passes,  a and  [3  approach  zero,  applying  heavy  smoothing  to  the  new  sam- 
ples. This  method  is  optimal  for  straight-line  tracks  but  must  be  modified  to  enable  the  filter 
to  follow  target  maneuvers. 


Maneuver-Following  Logic 

Benedict  and  Bordnert  noted  that  in  track-while-scan  systems  there  is  a conflicting  re- 
quirement between  good  noise  smoothing  (implying  small  o and  (3)  and  good  maneuver- 
following capability  (implying  large  o and  fi).  Although  some  compromise  is  always  required, 
the  smoothing  equations  should  be  constructed  to  give  the  "best"  compromise  for  a desired 


i 

| 


t 


i| 

i 


*A  L Qutglev.  IEEE  Internationa!  Conference  on  Radar  — Present  and  Future.  352*357,  1973. 
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noise  reduction.  Specifically,  since  the  variance-reduction  ratio  K,  defined  as  the  steady-state 
variance  in  the  filter  position  output  divided  by  the  variance  in  the  measured  position,  equals 


2a  2 + fi  (2  — 3<> ) 
a [4  - fi  - 2a  1 


(68) 


the  (a,  fi)  pair  should  be  chosen  to  satisfy  (68)  and  maximize  the  maneuver-following  capabili- 
ty. Benedict  and  Bordner*  defined  a measure  of  transient-following  capability  and  showed  that 
<«  and  (i  should  be  related  by 

«2 

(i  . (69) 

l — a 


Alternately  an  (<t , fi)  pair  satisfying  (69)  can  be  chosen  so  that  the  tracking  filter  will  follow-  a 
specified  g turn.  Cantrellt  developed  a method  of  determining  the  probability  that  the  target 
detection  will  fall  within  a correlation  region  centered  at  the  predicted  target  position  when  the 
target  is  performing  a specified  g-turn.  Then  the  («,  fi)  pair  yielding  the  smallest  correlation  re- 
gion should  be  used. 

The  trouble  with  the  preceding  method  is  that  if  high-^j  turns  are  followed  the  noise  per- 
formance is  rather  poor.  To  rectify  this  situation,  a turn  detector  employing  the  two  correlation 
regions  shown  in  Fig.  22  is  used.  If  the  detection  is  in  the  nonmaneuvering  correlation  region, 
the  filter  operates  as  usual,  a and  fi  being  reduced  according  to  (66)  and  (67).  Usually  it  is 
worthwhile  to  bound  a and  li  from  zero  by  assuming  a random  acceleration  Q(t)  ^ 0 
corresponding  to  approximately  a 1-g  maneuver.  When  the  target  falls  outside  the  inner  gate 
but  within  the  maneuver  gate,  a maneuver  is  declared  and  the  filter  bandwidth  is  increased  (<* 
and  fi  are  increased);  Quigley  and  Holmes*  increase  the  bandwidth  by  lowering  the  value  of  k 
in  (66)  and  (67).  To  avoid  the  problem  of  the  target  fading  and  a false  alarm  appearing  in  the 
large  maneuver  gate,  the  track  should  be  bifurcated  when  a maneuver  is  declared.  That  is,  two 
tracks  are  generated:  the  old  track  with  no  detection  and  a new  maneuvering  track  with  the  new 
detection  and  increased  bandwidth.  The  next  detection  is  used  to  resolve  the  ambiguity  and 
remove  one  of  the  tracks. 


MANEUVERING  GATE 


PREDICTED 
. POSITION 


Fig.  22  — Maneuver  and  nonmaneuver  gales  centered  ai 
ihe  target's  predicted  portion 


-A 


SMOOTHED 

POSITION 


L NON  MANEUVERING 
GATE 


*T.  R.  Benedict  and  G.  W.  Bordner,  I FEE  Trans  Automatic  Control  AC-7  iNu  41  37-32  (I9h2) 

tB.  II.  Cantrell.  "Behavior  of  » - li  Tracker  for  Maneuvering  targets  Under  Noise,  f alse  Target,  and  fade  Conditions." 
NRL  Report  7434.  Aug  1972. 

fA,  L.  Quigley  and  J.  E.  Holmes.  "The  development  of  Algorithms  for  the  Formation  and  Updating  of  Tracks,  Ad- 
mirably Surface  Weapons  Establishment.  WP-XBC-7512.  Portsmouth  P06  4A  A.  Nov  1975. 
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Cantrell  et  al 

.*  suggested  that  the  <»  — p filter  (described  by  (55),  (56),  and  (57))  be  made 

adaptive  by  adjucti 

->g  a and  [i  by 

_ . —2fc>n  T 

it  =1  — e 

(70) 

and 

13  =1  + e-2f"’° r ~2e~e'""1  cos((u07'a/i  - €2). 

' (71) 

in  which 

ojq  = 0 .5 1 ( k )l Pi  ( k ) | , 

(72) 

where 

px  (A)  =e  (A  -1 ) + (\  - e "’hT)  e (A)e  (A  —1  ) 

(73) 

and 

p2(A)  =e  T p2  (A  —1  ) + (1  -e  b’bT)  e(A)e(A) 

(74) 

with  £ being  the  damping  coefficient  (nominally  0.7),  T being  the  time  since  the  last  update, 
o>a  and  o>b  being  weighting  constants,  and  «(A)  being  the  error  between  the  measured  and 
predicted  positions  on  the  Ath  update.  The  basic  principle  of  the  filter  is  that  px  (A)  is  an  esti- 
mate of  the  covariance  of  successive  errors  and  (A)  is  an  estimate  of  the  error  variance. 
When  the  target  trajectory  is  a straight  line,  P\  (A)  approaches  0,  since  the  expected  value  of 
« (A)  is  0.  Thus  w0  approaches  0,  and  the  filter  performs  heavy  smoothing.  When  the  target 
turns,  (A)  grows,  since  the  error  e (A)  will  have  a bias  (either  positive  or  negative).  Thus  oi0 
grows,  and  the  filter  can  follow  the  target  maneuver. 

Another  solution  to  the  target-maneuvering  problem  is  due  to  Singert,  who  suggested  us- 
ing the  Kalman  filter  with  a realistic  target-maneuvering  model.  He  assumed  that  the  target 
was  moving  at  a constant  velocity  but  was  being  perturbated  by  a random  acceleration.  Since 
the  target  acceleration  is  correlated  in  time  (if  the  target  is  accelerating  at  time  /,  it  is  likely  to 
be  accelerating  at  time  t + Sr),  it  was  assumed  that  the  covariance  of  the  correlation  was 


/•(r)  = E\a(t)a(t+-))  -w*  (75) 

where  a(r)  is  the  target  acceleration  at  time  I,  <r}n  is  the  variance  of  the  target  acceleration,  and 
« is  the  reciprocal  of  the  maneuver  time  constant.  The  density  function  for  target  acceleration 
consists  of  delta  functions  at  ± /(max,  each  with  probability  7^.^,  a delta  function  at  0 with 
probability  Pn,  and  a uniform  density  between  —Am ax  and  .4max . For  this  density 


.2  _ 


(1 


+ 4Pmjx  -V- 


(76) 


lor  this  target  motion  Singert  then  calculated  the  state  transition  matrix  <!>U)  and  the  covari- 
ance matrix  Q<t).  thereby  specifying  the  Kalman-filter  solution.  He  generated  curves  which 
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give  the  steady-state  performance  of  the  filter  for  any  data  rate,  single-look  measurement  accu- 


racies, encounter  geometry,  and  class  of  maneuvering  targets. 

I 

Track  Initiation 

Detections  that  do  not  correlate  with  clutter  points  or  update  tracks  are  used  to  initiate 
new  tracks.  If  the  detection  does  not  contain  doppier  information,  the  new  detection  is  used  as 
the  predicted  position,  and  a large  correlation  region  must  be  used.  The  probability  of  false 
alarms  being  in  the  large  correlation  region  is  large;  hence  tracks  should  not  be  declared  firm 
until  a third  detection  (falling  within  a smaller  correlation  region)  is  obtained.  The  usual  initia- 
tion criteria  are  three  out  of  four  and  three  out  of  five.  The  possible  exceptions  are  when 
doppier  information  is  available  (so  that  a small  correlation  region  can  be  used  immediately)  or 
for  popup  (close)  targets  in  a military  situation. 

Quigley  and  Holmes*  suggest  using  a sequential  hypothesis-testing  scheme  for  initiating 
tracks.  When  a correlation  is  made  on  the  / th  scan,A(  is  added  to  the  likelihood;  and  when  a 
correlation  opportunity  is  missed,  A,-  is  subtracted  from  the  likelihood.  The  increment  A,  is  set 
by  the  state  of  the  tracking  system,  being  a function  of  the  closeness  of  the  association,  the 
number  of  false  alarms,  the  a-priori  probability  of  targets,  and  the  probability  of  detection. 
Although  this  method  will  inhibit  false  tracks  in  dense  detection  environments,  it  will  not 
necessarily  establish  the  correct  tracks.  The  proper  solution  will  probably  be  a method  of  gen- 
erating trial  tracks  using  detections  front  the  last  several  scans  and  then  eliminating  the  false 
tracks  in  an  easily  implemenable  manner. 

To  initiate  tracks  with  detections  front  unsynchronized  radars.  Trunk  et  al.t  suggest  that 
two  limes,  namely  T0  and  be  used.  If  time  T0  goes  by  between  track  updates,  the  tenta- 
tive track  should  be  dropped;  and  if  the  track  is  updated  after  a time  TF  after  initialization,  the 
tentative  track  should  be  made  firm.  Setting  TF  > TD  insures  that  three  detections  are  re- 
quired for  making  a track  firm. 

Firm  tracks  that  are  not  updated  in  30  or  40  seconds  are  usually  dropped. 

Correlation  Logic 

Several  procedures  will  now  be  given  for  associating  detections  with  tracks.  Of  special  in- 
terest are  the  conflict'ng  situations  of  multiple  tracks  competing  for  a single  detection  or  of 
multiple  detections  lying  within  a track’s  correlation  gate  (or  region). 

First  to  limit  the  number  of  detections  that  can  update  a track,  correlation  gates  are  used. 
A detection  can  never  update  a t rack  unless  it  lies  within  the  correlation  gate  which  is  centered 
tit  the  track's  predicted  position.  The  correlation  gate  should  be  defined  in  r n coordinates,  re- 
gardless of  what  coordinate  system  is  being  used  for  tracking.  Furthermore  the  gate  size 
should  be  a function  of  the  measurement  accuracy  K(i)  and  prediction  error  /’(/|r  —I  ) so  that 

' \ I Qini;le>  rml  I I Holmes.  The  development  of  Algorithms  lor  the  f ormation  and  Updating  ot  Tracks  " Ad- 
mirality  Surface  Weapons  I stahlishment.  Wl’-X  BC-75 1 2.  Portsmouth  POP  4 A V Nov.  I ')  7 s 

tCi  V Trunk,  J D Wilson.  H It  Cantrell.  J J.  Alter,  and  I I).  Queen.  "Moddlcalinns  U>  and  Preliminary  Results  for 
the  ADIT  System.'  NRl.  Report  8091.  Apr.  1077 
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the  probability  of  the  correct  detection  lying  with  the  gate  is  high  (at  least  0.99).  In  some 
tracking  systems,*  the  correlation  gate  is  fed  back  to  the  automatic  detector,  and  the  detection 
threshold  is  lowered  in  the  gate  to  increase  PD. 


When  several  detections  are  within  the  correlation  region,  the  usual  and  simplest  solution 
is  to  associate  the  closest  detection  with  the  track.  Specifically,  the  measure  of  closeness  is  the 
statistical  distance 


D 2 


(77) 


where  ( rp , Hp)  is  the  predicted  position,  (r,„.  »„,)  is  the  measured  position,  <t } is  the  variance 
of  /■„  - r„,  and  ( r } is  the  variance  of  Hn  - Since  the  prediction  variance  is  proportional 
to  the  measurement  variance,  arf  and  Of  are  sometimes  reptaced  by  the  measurement  vari- 
ances. Statistical  distance  rather  than  Euclidean  distance  must  be  used,  because  the  range  ac- 
curacy is  usually  much  better  than  the  azimuth  accuracy. 


Problems  associated  with  multiple  detections  and  tracks  are  illustrated  in  Fig.  23:  two 
detections  are  within  gate  1,  three  detections  are  within  gate  2,  and  one  detection  is  within  gate 
3.  Table  3 lists  all  detections  within  the  tracking  gate,  and  the  detections  are  entered  in  the 
order  of  their  statistical  distance  from  the  track.  Tentatively,  the  closest  detection  is  associated 
with  each  track,  and  then  the  tentative  associations  are  examined  to  remove  detections  which 
are  used  more  than  once.  Detection  8,  which  is  associated  with  tracks  l and  2,  is  paired  with 
the  closest  track  (track  1 in  this  case);  then  all  other  tracks  are  reexamined  to  eliminate  all  as- 
sociations with  detection  8.  Detection  7 is  tentative  associated  with  track  2;  a conflict  is  noted 
but  is  resolved  by  pairing  detection  7 with  track  2.  When  other  associations  with  detection  7 
are  eliminated,  track  3 has  no  associations  with  it  and  consequently  will  not  be  updated  on  this 
scan.  Thus  track  1 is  updated  by  detection  8,  track  2 is  updated  by  detection  7,  and  track  3 is 
not  updated. 


Fig.  23  ~ Example  of  ihe  problems 

caused  by  . multiple  detections  and 
tracks  in  close  vicinity 


An  alternate  strategy  is  to  always  pair  a detection  with  a track  if  there  is  only  one  correla- 
tion with  a track.  As  before,  ambiguities  are  removed  by  using  the  smallest  statistical  distance. 
Thus  track  3 in  the  example  is  updated  by  detection  7,  track  1 is  updated  by  detection  8,  and 
track  2 is  updated  by  detection  9. 


•x  R t ook.  Development  of  IADT  Tracking  Algorithm.  Johns  Hopkins  l ntverv.lv.  Vpp Physic  Laboratory. 
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Table  3 - Associaiion  Table  for  the  Example  shown  in  Pig.  23 


Track 

Number 

Closest  Association 

Second  Association 

Third  Association 

Detection 

Number 

D2 

Detection 

Number 

D2 

Detection  , 
Number 

D2 

l 

8 

1.2 

7 

4.2 

_ 

HI 

2 

8 

3.1 

7 

5.4 

9 

B 1 

3 

7 

6.3 

- 

- 

- 

mm 

Singer  and  Sea*  were  two  of  the  first  people  to  recognize  and  characterize  the  interaction 
between  the  correlation  and  track  update  functions.  Specifically,  three  distinct  situations  can 
occur:  the  track  is  not  updated,  the  track  is  updated  with  the  correct  return,  and  the  track  is 
updated  with  an  incorrect  return.  They  generalized  the  tracking  filter’s  error  covariance  equa- 
tions to  account  for  the  a priori  probability  of  incorrect  returns  being  correlated  with  the  track. 
This  permits  the  analytical  evaluation  of  tracking  accuracy  in  a multitarget  environment  which 
produces  false  correlations.  Furthermore,  using  the  generalized  tracking  error  covariance  equa- 
tion, they  optimized  the  filter  gain  matrix,  which  yielded  a new  minimum-error  tracking  filter 
for  multitarget  environments.  Also,  they  generated  a suboptimal  fixed-memory  version  of  this 
filter  to  reduce  computation  and  memory  requirements. 

A later  paper  by  Singer  el  al.f  uses  a-posteriori  correlation  statistics  based  on  all  reports  in 
the  vicinity  of  the  track.  Again  the  mathematical  structure  is  similar  to  the  Kalman  filter:  the 
state  equation  is  (58),  the  observation  equation  is  (59),  the  one-step  prediction  is  (60),  and  the 
corresponding  covariance  matrix  is  (61).  The  estimation  error  is  denoted  by  A" (/|  z ) = 
X(t)  — X(t\t')  and  has  mean  and  covariance  matrices  denoted  by  b(t\t')  and  The 

correlation  gate  size  and  shape  is  based  on  the  Mahalanobis  distance,*  and  it  is  assumed  nk 
sensor  reports  fall  within  the  gate  on  scan  A.  Included  in  the  number  nk  are  extraneous  re- 
ports whose  number  obeys  a Poisson  distribution  and  whose  positions  are  uniformly  distribut- 
ed within  the  gate.  The  smooth  estimate  is  given  by 

A'(fU)  = XU\i  -1 ) + AU).  (78) 

where  A (r)  is  chosen  to  minimize  the  noncentral  second  moment  of  the  filter  estimation  error. 

The  problem  is  solved  by  using  track  histories.  A track  history  a at  scan  A-  is  defined  by 
selecting,  for  each  scan  / < A,  a sensor  report  Y1'  (i)  w here  0 < y,  < nr  with  j\  =0 
corresponding  to  the  hypothesis  that  none  of  the  reports  belong  to  the  track.  The  number  of 
such  track  histories  is 

MA)  = n ('  + ",  >•  179) 


*K  A.  Singer  and  R.  G Sea,  li  t- K Trans  Automatic  Cnnirol  AC-18.  571-582  (1973). 

tR  A Singer,  R G.  Sea.  and  K B llouscwright,  IEF.K  Trans  Information  Theory  IT-20,  423-432  < 1974) 

>R  A Singer  and  A J Kan>uck.  Automation  7,  455-4h3  (1971) 
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Associated  with  each  history  <«  is  the  probability  p„  ( t ) that  the  history  u is  the  correct  one, 
given  observations  through  time  r (scan  A).  The  terms  6lt  (r|r-l)  and  Piy  (/!/-l)  are  the 
bias  and  covariance  of  the  estimation  error  X(r\t  —1  ),  given  observations  through  time  t — 1 
and  given  that  track  history  «'at  time  r — 1 is  the  (only)  correct  one.  Recursive  equations  are 
obtained  for  piy , bt< . and  Ptl ; then  it  is  shown  that  the  optimal  correction  vector  is  given  by 

Uk) 

A U)  = £ P„  U)bn  <'l'  “I  >•  (80) 

This  solution  not  only  minimizes  the  mean-square  error  but  also  is  an  unbiased  estimate. 

The  trouble  with  the  optimal  a-posteriori  filter  is  that  it  requires  a growing  memory. 
Hence  several  suboptimal  filters  were  suggested.  The  first  suboptimal  filter  considers  only  the 
last  /V  scans;  track  histories  which  are  identical  for  the  last  A' scans  are  merged.  The  second 
suboptimal  filter  only  considers  the  L nearest  neighbors  in  the  correlation  gate;  essentially  the 
gate  size  is  changed  to  limit  number  of  reports  to  L.  The  last  method  uses  both  techniques: 
considers  only  the  last  A' scans  and  restricts  number  of  reports  on  any  scan  to  L. 

Simulations  were  run  to  compare  the  optimal  and  suboptimal  a-posteriori  filters,  the  op- 
timal and  suboptimal  a-priori  filters,  and  the  Kalman  filter.  Some  of  the  results  are  summarized 
in  Figs.  24  and  25.  In  Fig.  24  the  filter  variance  normalized  by  the  theoretical  (perfect- 
correlation)  Kalman-filter  variance  is  plotted  for  several  filters.  As  a class  the  a-posteriori  filters 
provide  better  performance  than  the  other  filters.  However,  for  high  density  of  false  reports 
(4 =0.1),  the  a-posteriori  filter  is  30  times  worse  than  predicted  by  the  standard 
Kalman-filter  approach.  Thus  the  standard  approach  should  never  be  used  in  dense-target  (or 
false-target)  environments.  Figure  25  gives  the  probability  of  making  a false  correlation. 
Again  the  a-posteriori  filters  provide  the  best  performance. 

Stein  and  Blackman*  have  proposed  a maximum-likelihood  approach  similar  to  Sittlert  for 
solving  the  multitarget  correlation  problem.  Their  approach  is  unified  in  that  they  consider  the 
total  correlation-track  problem  which  includes  track  initiation,  confirmation,  gating,  and  dele- 
tion logic.  They  compare  their  results  with  a standard  approach  and  show  significant  perfor- 
mance improvements.  However  this  author  wonders  how  complicated  the  method  is  to  imple- 
ment and  what  their  improvement  is  relative  to  a more  sophisticated  approach  such  as  that  of 
Singer  et  al.t. 

Radar  Integration 

There  are  many  ways  of  integrating  (combining)  radar  detections  from  multiple  radars 
into  a single  system  track  file.  The  type  of  radar  integration  that  should  be  used  is  a function 
of  the  radar’s  performance  and  its  environment.  Although  no  firm  rules  can  be  generated, 
several  methods  and  some  general  rules  are  as  follows: 

• Track  selection.  Generate  a track  with  each  radar,  and  choose  one  of  the  tracks  for 
the  system  track.  The  only  advantage  of  this  method  is  that  it  is  the  simplest 
method  to  implement. 


*J  J Stein  ami  S S.  Blackman.  IEEE  Trans.  Aerospace  and  Electronic  Systems  AES-11.  1207-1217  (1975) 
♦R  W Sittler.  II* I E Trans  Military  Electronics  ME-8,  125-139  (1964) 

R \ Singer.  R (»  Sea.  and  K B House* right.  IEEE  Tranv  Information  Theorv  IT-20.  423-432  <14'4» 
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• Average  track.  Generate  a track  with  each  radar,  and  weight  the  tracks  to  form  a 

system  track.  The  method  can  be  applied  when  man>  radars  are  providing  unsyn- 
chronized radar  date. 

• Augmented  track.  Generate  a track  with  each  radar,  choose  one  as  the  system  track, 

but  also  use  selected  detections  from  other  radars  to  update  the  system  track.  This 
method  should  be  used  when  one  radar  provides  substantially  better  data  than  oth- 
er radars.  Detections  from  other  radars  should  be  used  when  the  primary  radar 
misses  some  detections  or  when  a target  maneuver  is  declared. 

• Average  detection.  Average  all  detections,  and  use  the  average  to  form  a system 

track.  This  method  should  be  used  when  many  radars  are  providing  detections 
essentially  at  the  same  instant. 

• Merged  detections.  Use  all  detections  to  update  the  system  track;  tracks  may  or  may 

not  be  initiated  using  all  detections.  Theoretically  this  method  provides  the  most 
information;  that  is,  if  the  detections  are  properly  weighted,  this  method  always 
provides  the  best  performance.  However  care  must  be  taken  so  that  bad  data  do 
not  corrupt  good  data.  Thus  this  method  should  be  used  when  the  radars  are  sup- 
plying data  of  comparable  accuracy. 

The  most  important  advantage  of  radar  integration  is  provide  the  tracking  information  in 
one  central  source.  Radar  integration  also  provides  improved  track  continuity  and  improved 
tracking  performance  on  maneuvering  targets.  Little  improvement  is  obtained  in  track- 
initiation  times,  since  in  practice  almost  always  one  radar  will  detect  and  establish  a track  be- 
fore any  other  radar  can  provide  some  detections. 

Either  the  <*  — /i  fillers  or  the  Kalman  filler  can  be  used  when  the  radars  are  at  the  same 
location.  However  Trunk  and  Wilson*  indicate  that  the  Kalman  filter  must  be  used  to  provide 
triangulation  effects  when  the  radars  are  in  different  locations.  Various  methods  for  multiple- 
site  correlation  are  also  discussed  by  Singer  and  Kanyuckt- 


CONCLUDING  REMARKS 

The  problems  involved  in  coherent  processing  have  received  the  greatest  attention  in  this 
survey.  Presently  the  trend  appears  to  be  toward  a digital  implementation  of  the  adaptive  pro- 
cessing algorithms:  direct  open-loop  calculation  of  canceler  weights  and  numerical  inversion  of 
the  covariance  matrix  for  adaptive  arrays.  It  can  be  expected  that  such  systems  will  be  built 
and  tested  during  the  next  several  years. 

The  area  of  noncoherent  processing  has  been  studied  intensively  since  1940.  The  em- 
phasis in  later  years  has  been  on  techniques  to  limit  the  number  of  false  alarms  (while 
suffering  only  a small  target-sensitivity  loss)  so  as  not  to  overload  the  tracking  system.  Many 
systems  have  been  built  and  tested,  and  not  much  more  research  seems  necessary  in  this  area. 


'Ci  V Trunk  and  J.  I)  Wilson.  "Tracking  fillers  lor  Multiple-Platform  Radar  Imegralion."  NRL  Report  8087.  Dec 
1976, 

t!<  A Singer  and  A J Kun>uck.  Automatic.!  7.  457  46.1  (1971) 
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During  the  last  several  years  much  progress  has  been  made  in  track-while-scan  systems. 
This  work  has  given  guidance  on  the  important  problem  of  track-detection  correlation  in  a 
dense  multitarget  environment.  The  major  problem  still  needing  a solution  is  that  of  track  ini- 
tiation in  a dense  environment. 

One  problem,  which  has  received  little  attention  so  far  but  which  will  receive  more  atten- 
tion, is  that  of  adaptively  controlling  the  surveillance  radar.  Problems  of  interest  are:  when 
should  frequency  and/or  polarization  diversity  be  used,  when  and  where  should  various  radar 
modes  be  used,  and  how  should  the  signal  processing  be  reconfigured  to  cope  with  a changing 
environment?  Future  radars  will  have  more  flexibility,  and  their  control  will  become  extreme- 
ly important. 
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